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Universal quantum computation on decoherence-free subspaces and subsystems (DFSs) is ex- 
amined with particular emphasis on using only physically relevant interactions. A necessary and 
sufficient condition for the existence of decoherence-free (noiseless) subsystems in the Markovian 
regime is derived here for the first time. A stabilizer formalism for DFSs is then developed which 
■ allows for the explicit understanding of these in their dual role as quantum error correcting codes. 

Conditions for the existence of Hamiltonians whose induced evolution always preserves a DFS are 
derived within this stabilizer formalism. Two possible collective decoherence mechanisms arising 
from permutation symmetries of the system-bath coupling are examined within this framework. It 
5_l , is shown that in both cases universal quantum computation which always preserves the DFS {natural 

Oh fault-tolerant computation) can be performed using only two-body interactions. This is in marked 

' contrast to standard error correcting codes, where all known constructions using one or two-body 

interactions must leave the codespace during the on-time of the fault-tolerant gates. A further 
consequence of our universality construction is that a single exchange Hamiltonian can be used to 
perform universal quantum computation on an encoded space whose asymptotic coding efficiency 
is unity. The exchange Hamiltonian, which is naturally present in many quantum systems, is thus 
' asymptotically universal. 

\o ■ 

O ' I. INTRODUCTION 

-xf : 

The discovery that information encoded over quantum systems can exhibit strange and wonderful computational 
jl],^| and information theoretic j^^] properties has led to an explosion of interest in understanding and exploiting the 
"quantumness" of nature. For the use of quantum information to progress beyond mere theoretical constructs into 
the realm of testable and useful implementations and experiments, it is essential to develop techniques for preserving 
quantum coherences. In particular, the coupling of a quantum system to its environment leads to a process known 
as decoherence, in which encoded quantum information is lost to the environment. In order to remedy this problem 
active quantum error correction codes (QECCs) |3|-|9| have been developed, by analogy with classical error correction. 
These codes encode quantum information over an entangled set of codewords, the structure of which serves to preserve 
the quantum information, when used in conjunction with a frequently recurring error correcting procedure. It has 
been shown that when the rate of decoherence is below a certain threshold, fault tolerant quantum information 
manipulation is possible fllo|-[l2f . Since it is believed that there are no systems for which the decoherence mechanism 
entirely vanishes, QECCs will be essential if quantum information manipulation is to become practical. 

An alternative approach has been proposed and developed recently, in which the central motivation is the desire to 
reduce the effect of a specific decoherence mechanism. This is the decoherence-free subspace (DFS) approach (also 
referred to as "error avoiding", or "noiseless" quantum codes) [^3|^26|. In contrast to the active mode of QECC, 
DFS theory can be viewed as providing a passive approach, where a specific symmetry of the system-bath coupling 
is employed in order to seek out a quiet corner of the system's Hilbert space which does not experience decoherence. 
Information encoded here over a subspace of (usually entangled) system states is robust against a specific form of 
decoherence. We shall refer to this as the "DFS supporting decoherence mechanism". When this is the dominant 
form of decoherence in the physical system, there are major gains to be had by operating in the DFS. Previous work 
has shown that collective decoherence of the type experienced in condensed phase systems at low temperatures can 
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be successfully eliminated in this way 27 28j] . Further research showed that DFSs are robust to perturbing error 
processes [^l],p4|, and are thus ideally suited for concatenation in a QECC p2| . 

The motivating goal behind the DFS approach is to use symmetry first. Thus, one first identifies a DFS for the major 
sources of decoherence, via the symmetry of the interaction with the environment. One then proceeds to use the DFS 
states as a basis for a QECC which can deal with additional perturbing error processes. In order for this scheme to be 
credible, DFSs must support the ability to perform universal quantum computation on the encoded states. Towards 



this end, certain existential results 29 have been derived showing that in principle universal quantum computation 
can be performed on any DFS. Constructive results for a set of universal quantum gates on a particular class DFSs 
were subsequently constructed in |23] using known QECC constructions. However, these gates were constructed in 



1 



such a way that during the operation of the gate, states within a DFS are taken outside of this subspace. Thus these 
gates would necessarily need to operate on a timescale faster than the DFS supporting decoherence mechanism, in 
order to be applied efficiently to a concatenated DFS-QECC scheme. R Similarly, a universal computation result on 
DFSs for atoms in cavities was recently presented by Beige et al. in [ pOpifl . It assumes that the interaction driving a 
system out of the DFS is much weaker than the coupling of non DF-states to the environment. It is then possible to 
make use of an environment-induced quantum Zeno effect. In order to make use of the robustness condition without 
resorting to gates which can be made faster than the main DFS supporting decoherence mechanism, one would prefer 
to explicitly construct a set of Hamiltonians which can be used to perform universal quantum computation, but 
which never allow states in the DFS to leak out of the DFS. Imperfections in these gates may be dealt with by the 
concatenation technique of J22| (see also |2j|). 

In addition, one would, from a practical standpoint, like to use Hamiltonians which involve at most two-body 
interactions (under the assumption that any three-body interactions will be weak and not useful for operations 
which must compete with the decoherence rate). In p2| such Hamiltonians were used for the important decoherence 
mechanism known as "collective decoherence" , on a system of 4 physical qubits. In collective decoherence the bath 
cannot distinguish between individual system qubits, and thus couples in a collective manner to the qubits. The 
corresponding two-body Hamiltonians used to implement universal quantum computation are those that preserve the 
collective symmetry: the exchange interaction between pairs of qubits. The first and main purpose of this paper is 
therefore to extend the constructive results obtained in |3^] to other forms of collective decoherence and to larger DFSs. 
Two different forms of collective decoherence are considered here, and constructive results are obtained for these on 
DFSs of arbitrary numbers of qubits. These results have implications that extend far beyond the problem of dealing 
with collective decoherence. Since they imply that the exchange interaction by itself is sufficient to implement universal 
quantum computation on a subspace, it follows that using encoded (rather than physical) qubits can be advantageous 
when resources for physical operations are limited. After all, the standard results for universal quantum computation 
employ either arbitrary single-qubit operations in addition to a non-trivial two-qubit gate (e.g., a controlled-NOT), or 
at least two non-commuting two-qubit Hamiltonians []33|-[37|]. These issues will be explored in a separate publication. 

Previous work established that DFSs correspond to the degenerate component of a QECC [^,^8|. A second 
purpose of this work is to present new results on a recently discovered generalization of DFSs, which has been 
termed "noiseless subsystems", and arises from a theory of QECC for general decoherence mechanisms |3^,^(J. In 
line with our previously established terminology pl[ we will refer to these as "decoherence-free subsystems" , where 
we take the term "decoherence" to mean both dephasing (T 2 ) and dissipation (Ti). Essentially, the generalization 
corresponds to allowing for information to be encoded into states transforming according to ar6i£rar?/-dimensional 
irreducible representations (irreps) of the decoherence-operators' algebra, instead of just one-dimensional irreps as in 
the decoherence-free subspace case (we will present precise definitions later in this paper). These results all arise from 
a basic theorem on algebras that are closed under the Hermitian conjugation operation ("f -closed algebras"), and 
thereby unify the role of symmetry in both decoherence-free subspaces and quantum error correction. In this paper 
we extend the decoherence-free subsystem concept to situations governed by essentially non-f-closed evolution. Such 
situations arise from non-Hermitian terms in the system-bath interaction, which may occur, e.g., in generalized master 
equation and conditional Hamiltonian representations of open quantum dynamics J4~i] . In particular, we derive an if 
and only if (iff) condition for the existence of decoherence-free subsystems with dynamics governed by a semigroup 
master equation. This is important because it is well known in decoherence-free subspace theory that such non-f- 
closed evolution can support different DFSs than in the f-closed case. A similar result is now shown here to hold for 
the decoherence-free subsystems. 

Existential results for universal quantum computation on decoherence-free subsystems also exist [ [i"2"f . The universal 
quantum computation results we obtain in this paper extend beyond decoherence-free subspaces: we show how to 
achieve constructive universal quantum computation on the decoherence-free subsystems supported under collective 
decoherence. This most significant achievement of our paper settles the question of universal quantum computation 
under collective decoherence using realistic Hamiltonians. 

Another aim of this paper is to elucidate the close link between DFS and QECC. In ^||,[58) it was shown that DFSs 
are in fact maximally degenerate QECCs. This result was derived from the general condition for a code to be a QECC 
H . A very fruitful approach towards QECC has been the stabilizer formalism developed in || which led to the theory 
of universal fault-tolerant computation on QECCs (43). In |^3| we considered DFSs as abelian stabilizer codes. Here 
we generalize the stabilizer-framework to non-abelian stabilizers, and show that in general DFSs are stabilizer-codes 
that protect against errors in the stabilizer itself. This perspective allows in return to view QECCs as DFSs against 



1 Note that QECC fault-tolerant gates are also required to operate faster than the decoherence time of the main error process. 
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a certain kind of errors, and establishes a kind of duality of QECCs and DFSs. 

The paper is structured as follows: In Section || we review decoherence-free subsystems and place them into the 
context of the Markovian master equation. For decoherence-free subspaces this has been done in ||l],^4| . These earlier 
results are therefore generalized here to subsystems. In Section ^ we introduce a generalized stabilizer-formalism 
for DFS, and connect to the theory of stabilizers on QECC developed in |^]. This allows us to treat DFS and QECC 
within the same framework. It also sheds some light on the duality between DFS and QECC, in particular on the 
performance of a DFS viewed as a QECC and vice versa. In Section [V we deal with universal computation on DFS 
within both the stabilizer-framework and the representation-theoretic approach. We derive fault-tolerance properties 
of the universal operations. In particular, we show how to obtain operations that keep the states within a DFS during 
the entire switching-time of a gate. Further we define the allowed compositions of operations and review results on 
the length of gate sequences in terms of the desired accuracy of the target gates. In Section [V| we introduce the model 
of collective decoherence. Section VI explicitly deals with the abelian case of weak collective decoherence in which 
system-bath interaction coupling involves only a single system operator. Stabilizer and error-correcting properties 
are developed for this case, and it is shown how universal computation ca n be achieved. The same is done for the 
non-abelian and more general case of strong collective decoherence in Section VTJ . For both weak and strong collective 
decoher ence we show how to fault-tolerantly encode into and read out of the respective DFSs. Finally, we analyze in 
Section VIII how to concatenate DFSs and QECCs to make them more robust against perturbing errors (as proposed 
in and show how the universality results can be applied to achieve fault- tolerant universal computation on these 
powerful concatenated codes. We conclude in Section IX. Derivations and proofs of a more technical nature are 
presented in the Appendix. 



II. OVERVIEW OF DECOHERENCE-FREE SUBSPACES AND SUBSYSTEMS 



A. Decoherence-Free Subspaces 



Consider the dynamics of a system S (the quantum computer) coupled to a bath B via the Hamiltonian 

H = H 5 (8)I s + Is«)H B + H / , (2.1) 

where H5 (H^) [the system (bath) Hamiltonian] acts on the system (bath) Hilbert space Hs (T~(-b), Is (Is) is the 
identity operator on the system (bath) Hilbert space, and H/, which acts on both the system and bath Hilbert spaces 
Ti.s ®7~Ib, is the interaction Hamiltonian containing all the nontrivial couplings between system and bath. In general 
H/ can be written as a sum of operators which act separately on the system (S a 's) and on the bath (B a 's): 

H / = ^S Q ®B a . (2.2) 



In the absence of an interaction Hamiltonian (H/ = 0), the evolution of the system and the bath are separately 
unitary: U(i) = cxp[— iHt] = cxp[— iHst] (g> exp[— iRst] (we set U = 1 throughout). Information that has been 
encoded (mapped) into states of the system Hilbert space remains encoded in the system Hilbert space if H/ = 0. 
However in the case when the interaction Hamiltonian contains nontrivial couplings between the system and the bath, 
information that has been encoded over the system Hilbert space does not remain encoded over solely the system 
Hilbert space but spreads out instead into the combined system and bath Hilbert space as the time evolution proceeds. 
Such leakage of quantum information from the system to the bath is the origin of the decoherence process in quantum 
mechanics. 

Let Tls be a subspace of the system Hilbert space with a basis \i). The evolution of such a subspace will be unitary 
@H if and only if (i) 

S a \i) =c a \i), c a e€ (2.3) 

for all \i) £ H.s and for all S a , (ii) Us does not mix states within the subspace with states that are outside of the 
subspace ((j'|Hs|?) = for all \i) in the subspace and all \j'} outside of the subspace: H5 = H5 H^) and (iii) 
system and bath are initially decoupled p(0) = Ps(0) Pb(0)- We call a subspace of the system's Hilbert space which 
fulfills these requirements a decoherence-free subspace (DFS). 

The above formulation of DFSs in terms of a larger closed system is exact. It is extremely useful for finding 
DFSs, providing often the most direct route via simple examination of the system components of the interaction 
Hamiltonian. In practical situations, however, the closed-system formulation of DFSs is often too strict. This is 
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because the closed-system formulation incorporates the possibility that information which is put into the bath will 
back-react on the system and cause a recurrence. Such interactions will always occur in the closed-system formulation 
(due to the the Hamiltonian being Hermitian). However, in many practical situations the likelihood of such an event 
is extremely small. Thus, for example, an excited atom which is in a "cold" bath will radiate a photon and decohere 
but the bath will not in turn excite the atom back to its excited state, except via the (extremely long) recurrence time 
of the emission process. In these situations a more appropriate way to describe the evolution of the system is via a 
quantum dynamical semigroup master equation |4~I||4E| ] . By assuming that (i) the evolution of system density matrix 
is a one-parameter semigroup, (ii) the system density matrix retains the properties of a density matrix including 
"complete positivity" , and (iii) the system and bath density matrices are initially decoupled, Lindblad |Q has shown 
that the most general evolution of the system density matrix ps(t) is governed by the master equation 

1 M 

LflM*)] = 2 ^ aafJ ([ F -'Ps(t)F^] + [F a p s (t),Fj]) (2.4) 

a, 13=1 

where Hg is the system Hamiltonian, the operators F Q constitute a basis for the M-dimensional space of all bounded 
operators acting on Tts, and a Q/ 3 are the elements of a positive semi-definite Hermitian matrix. As above, let lis be 
a subspace of the system Hilbert space TLs with a basis The evolution over such a subspace is then unitary [plj iff 

F a |i) = Ca|t>, c Q GC (2.5) 



for all \i) and for all F Q . While this condition appears to be identical to Eq (2.3), there is an important difference 
between the S Q 's and the F a 's which makes these two decoherence-freeness conditions different. In the Hamiltonian 
formulation of DFSs, the Hamiltonian is Hermitian. Thus the expansion for the interaction Hamiltonian Eq. (^^) can 
always be written such that the S Q are also Hermitian. On the other hand, the F Q 's in the master equation, Eq. (|2.4|) , 
need onl y be bounded operators acting on Hs and thus the F a 's need not be Hermitian. Because of this difference, 



Eq. (2.5) allows for a broader range of subspaces than Eq. (^3J). For example consider the situation where there are 
only two nonzero terms in a master equation for a two-level system, corresponding to Fi = er_ and F2 = a z where 
er_ = |0)(1| and a z — |0)(0| — |1)(1| (e.g., cooling with phase damping). In this case there is a DFS corresponding to 
the single state |0). In the Hamiltonian formulation, inclusion of Si = er_ in the interaction Hamiltonian expansion 
Eq. (|2.2|) would necessitate a secon d te rm in the Hamiltonian with S2 = <r_, along with the S z — a z as above. For 



this set of operators, however, Eq. (2.3) allows for no DFS 



B. Decoherence-Free Subsystems 

If one desires to encode quantum information over a subspace and requires that this information remains 



decoherence-free, then Eqs. (2.3),(p.q) provide necessary and sufficient conditions for the existence of such DFSs. 
The notion of a subspace which remains decoherence-free throughout the evolution of a system is not, however, 
the most general method for providing decoherence-free encoding of information in a quantum system. Recently, 
Knill, Laflamme, and Viola B9] have discovered a method for decoherence-free coding into subsystems instead of into 
subspaces. 

Decoherence-free subsystems p9M3 , ffl are most easily presented in the Hamiltonian formulation of decoherence. 
Let A denote the associative algebra formed by the system Hamiltonian H5 and the system components of the 
interaction Hamiltonian, the S Q 's. To simplify our discussion we will assume that the system Hamiltonian vanishes. 
(It is easy to incorporate the system Hamiltonian into the S Q 's when one desires that the system evolution preserves 
the decoherence-free subsystem.) We also assume that the identity operator is included as So = Is and B = 1b- 
This will have no observable consequence but allows for the use of an important representation theorem. A consists 
of linear combinations of products of the S a 's. Because the Hamiltonian is Hermitian the S Q 's must be closed under 
Hermitian conjugation: A is a f-closed operator algebra. A basic theorem of such operator algebras which include the 
identity operator states that, in general, A will be a reducible subalgebra of the full algebra of operators on Hs [(47). 
This means that the algebra is isomorphic to a direct sum of dj x dj complex matrix algebras, each with multiplicity 
nr. 

A*i(§)I nj ®M(dj,C). (2.6) 

JGJ 
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Here J is a finite set labeling the irreducible components of A, and A4(dj,€) denotes a dj x dj complex matrix 
algebra. It is also useful at this point to introduce the commutant A' of A. This is the set of operators which commutes 
with the algebra A, A' — {X : [X, A] = 0, VA e A}. They also form a f-closed algebra, which is reducible to 



A' = ®M(nj,€)®I dj 

JEJ 

over the same basis as A in Eq. (|2.6| ). 

The structure implied by Eq. (|2.6|) is illustrated schematically as follows, for some system operator S c 



(2.7) 



J = 1 



J = 2 



(2. 



where a typical block with given J has the structure: 



J given: 









M a 







d, - 1 



p! : • • ■ dj-1 



M a 



A = 



A = 1 



A = 71.7 - 1 



(2.9) 



Associated with this decomposition of the algebra A is the decomposition over the system Hilbert space: 



(2.10) 



Decoherence -free subsystems are defined as the situation in which information is encoded in a single subsystem s pace 
<D nj of Eq. ( |2.10|) (thus the dimension of the decoherence-free subsystem is n,j). The decomposition in Eq. ( |2.6| ) 
reveals that information encoded in such a subsystem will always be affected as identity on the subsystem space C™ J , 
and thus this information will not decohere. It should be noted that the tensor product nature which gives rise to 



the name subsystem in Eq. (2.6) is a tensor product over a direct sum, and therefore will not in general correspond 
to the natural tensor product of qubits. Further, it should be noted that the subsystem nature of the decoherence 
implies that the information should be encoded in a separable way. Over the tensor structure of Eq. (2.10) the density 
matrix should split into two valid density matrices: ps(0) — p ® 7 where p is the decoherence-free subsystem and 7 
is the corresponding component of the density matrix which does decohere. Finally it should be pointed out that not 
all of the subsystems in the different irreducible representations can be simultaneously used: (phase) decoherence will 
occur between the different irreducible components of the Hilbert space labeled by J <G J . For this reason, from now 
on we restrict our attention to the subspace defined by a given J. 

Decoherence-free subspaces are now easily connected to decoherence-free subsystems. Decoherence-free subspaces 
correspond to decoherence-free subsystems possessing one- dimensional irreducible matrix algebras: M(1,<C). The 
multiplicity of these one-dimensional irreducible algebras is the dimension of the decoherence-free subspaces. In fact 
it is easy to see how the decoherence-free subsystems arise out of a non-commuting generalization of the decoherence- 
free subspace conditions. Let {lA^)}, 1 < A < nj and 1 < fj, < dj, denote a subspace of TLs with given J. Then the 
condition for the existence of an irreducible decomposition as in Eq. (2.6) is 



dj 

E 



(2.11) 



for all S Q , A and p. Notice that M^^i jQ is not dependent on A, in the same way that c a in Eq. (2.3) is not the same 
for all \i) (there p = 1 and fixed). Thus for a fixed A, the subspace spanned by |A M ) is acted upon in some nontrivial 
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way. However, because M MM / !Q , is not dependent on A, each subspace denned by a fixed fi and running over A is acted 
upon in an identical manner by the decoherence process. 

At this point it should be noted that the generalization of the Lindblad master equation Eq. ( |2.4| ) with a decoherence- 
free subspace to the correspond ing master equation for a decoherence-free system is not trivial. This is because, as 
above, the F a operators in E q. (Bj ) are (for all practical purposes) not required to be closed under conjugation. The 
representation theo rem Eq. ( |2 .6|) is hence not directly applicable. We will show, however, that the master equation 
analog of Eq. (2.11) 

dj 

F Q |A M ) = J2 M^, a \X^) (2.12) 

provides a necessary and sufficient condition for the preservation of decoherence-free subsystems. 

As above, we consider a subspace of the system Hilbert space spanned by |A M ), with 1 < A < nj and 1 < /x < dj. 
Our notation will be significantly simpler if we explicitly write out the formal tensor product over this subspace: 
|A M ) = | A) ® \fi). In the subsystem notation, we claim that the decoherence-free subsystem condition is 

F a |A)<g>|/i) = |A)®M a |/i). (2.13) 

A proper decomposition of the system Hilbert space requires, as noted above, that the system density matrix is a 
tensor product of two valid (Hermitian, positive) density matrices: 

Ps(0)= £ m'(Q)7pM'(0)|A M )<A^|=p(Q)® 7 (0), (2-14) 

AA' ,fifi f 

where p(0) contains the information which will remain decoherence-free, and 7(0) is an arbitrary but valid density 
matrix. 

In general the operators F a will not be decomposable as a single tensor product corresponding to p(0) ® 7(0). 
Rather, they will be a sum over such tensor products, cor responding to an expansion over an operator basis: F a = 
J2 P ® Mg. The decohering generator of evolution (2.4) thus becomes 

1<d\ps(P)] = \ E <M E ( 2N ^(°) N ? ® M£ 7 (0)Mj - N«t N >(0) ® M*M£ 7 (0) 

a/3 pq 

- p(0)N^NP ® 7 (0)M^M^ . (2.15) 



Tracing over the 7 component, and using the cyclic nature of the trace allows one to factor out a common 
Tr 7 (Mg 7 (0)M| t ), yielding: 



pq 



T^{L D [p s (Q)} = \ E aa0<£ (2NP P (0)N^ - N^N^(O) - p(0)NjNS 



The evolution of the p component of the density matrix thus satisfies the standard master equation (2.4), for which 
it is known that the evolution is decoherence-free [^lj if and only if 

N« |A) = c Q>g |A) Va. (2.16) 

This implies that the necessary and sufficient condition for a decoherence-free subsystem is 

f q = e C",? 1 ® M i = 1 ® E c °.9 M " = 1 M «' ( 2 - 17 ) 



which is the claimed generalization of the Hamiltonian condition of decoherence-free subsystems, Eq. (2.12). 

We will use the acronym DFS to denote both decoherence-free subsystems and their restriction, decoherence-free 
subspaces, whenever no confusion can arise. When we refer to DF subspaces we will be specifically referring to the 
one-dimensional version of the DF subsystems. 
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III. THE STABILIZER FORMALISM AND ERROR-CORRECTION 



In the theory of quantum error correcting codes (QECCs) it proved fruitful to study properties of a code by 
considering its stabilizer S. This is the group formed by those system operators which leave the codewords unchanged, 
i.e., they "stabilize" the code. Properties of stabilizer codes and the theory of quantum computation on these stabilizer 
codes have been developed in [^3|. In the framework of QECCs, the stabilizer allows on the one hand to identify the 
errors the code can detect and correct. On the other hand it also permits one to find a set of universal, fault-tolerant 
gates by analyzing the centralizer of S, defined as the set of operations that commute with all elements in S (equal 
to the normalizer - the set of operations that preserve S under conjugation - in the case of the Pauli group). In the 
context of QECCs, the stabilizer S is restricted to elements in the Pauli-group, i.e., the group of tensor products of 
I,X, Y, Z, and is a finite abelian group. 

The extension of stabilizer theory yields much insight into DFSs. We do this here by defining a non-abelian, and 
in certain cases infinite stabilizer group. The observation that DFSs are highly degenerate QECCs will appear 
naturally from this formalism. Such a generalized stabilizer has already been defined in previous work dealing with 
decoherence-free subspaces j}^], and its normalizer shown there to lead to identification of local gates for universal 
computation. A key consequence of this approach was the observation that the resulting gates do not take the system 
out of the DFS during the entire switching time of the gate. 

We now review and extend the results in j3^| to analyze the error detection and correction properties of DFSs and 
QECCs. We shall incorporate DFSs and QECCs into a unified framework, similarly to the representation-theoretic 
approach of J3^,^]. The question of performing quantum computation on a specific DFS will be addressed in the 
next section. 



A. The Stabilizer - General Theory 



An operator S is said to stabilize a code C if 

|#)eC iff S|*) = |*> VSeS. (3.1) 

The set of operators {S} form a group S, known as the stabilizer of the code |43|. Clearly, S is closed under 
multiplication. In the theory of QECC the stabilizers that have been studied are subgroups of the Pauli-group (tensor 
products of I, X,Y,Z ). Since any two elements of the Pauli group either commute or anticommute, a subgroup 
(and in particular the stabilizer), in this case is always abelian P3| . The code is thus the common eigenspace of the 
stabilizer elements with eigenvalue 1. 

In general an error- process can be described by the Kraus operator-sum formalism |4^,|^] : p — > A^pAj, . The 
Kraus-operators A p can be expanded in a basis {E a } of "errors". The standard QECC error- model assumes that 
errors affect single qubits independently. Therefore the theory of QECCs has focused on searching for codes that 
make quantum information robust against 1, 2,... or more erroneous qubits. Detection and correction procedures must 
then be implemented at a rate higher than the intrinsic error rate. In the QECC error-model, the independent errors 
are spanned by single-qubit elements (X, Y, Z). An analysis of the error-correction properties can then be restricted 
to correction of combinations of these basic errors (which are also members of the Pauli group) acting on a certain 
number of qubits simultaneously. 

The distance d of a QECC is the number of single-qubit errors that have to occur in order to transform one 
codeword in C to another codeword in C. An error E is detectable if it takes a codeword to a subspace of the 
Hilbcrt space that is orthogonal to the space spanned by C (this can be observed by a non-perturbing orthogonal von 
Neuman measurement). A distance d code can detect up to d — 1 errors. In order to be able to correct an error on 
a certain codeword the error (up to a degenerate action of different errors) also needs to be identified, so that it can 
be undone. Hence errors on different codewords have to take the codewords to different orthogonal subspaces. The 
above translates to the QECC-condition S: 

A QECC C can correct errors £ — {E Q } if and only if 

= CapSij VEajEfi e £. (3.2) 

The stabilizer of a QECC allows identification of the errors which the code can detect and correct M. Two types 
of errors can be dealt with by stabilizer codes: (i) errors E^E^ ^ I that anticommute with an S € S, and (ii) 
errors that are part of the stabilizer (E Q e S ). It is straightforward to see that both (i) and (ii) imply the QECC 
condition Eq. (|J). For case (i), if eJe 3 -S = -SEjEj, then (^lE^E^) = (* i |E] 3 E a S|1' i ) = -(%|SE^E a |$i) = 
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-(tfj-IE^Ealtfi). Hence (*.,-|E£E a |¥i) = and Cap — 5 a p. Errors of type (ii), E Q G S, leave the codewords unchanged 
and therefore trivially lead to Eq. ( ft.2| ). The first class, (i), are errors that require active correction. The second 
class, (ii), are "degenerate" errors that do not affect the code at all. QECCs can be regarded as (passive) DFSs for 
the errors in their stabilizer [^3|. Conversely, being passive, highly degenerate codes J22|, DFSs can be viewed as a 
class of stabilizer codes that protect against type (ii) errors [i.e., where the are linear combinations of elements 
generated (under multiplication) by the stabilizer], and against the (usually small) set of errors that anticommutc 
with the DFS stabilizer. The stabilizer thus provides a unified tool to identify the errors that a given code can deal 
with, as a DFS and as a QECC. An analysis of the properties of DFSs with a stabilizer in the Pauli-group has been 
carried out in |25(] . 



B. DFS-Stabilizer 



Most of the DFSs stemming from physical error-models will not have a stabilizer in the Pauli-group, i.e., they are 
non-additive codes. The stabilizer may even be infinite. In particular, the codes obtained from a noise model where 
errors arise from a symmetric coupling of the system to the bath and that form the focus of this paper, are of this 
type. 

As discussed in the previous section, a DFS is completely specified by the condition: 



S q |m) ® | A) - |m) ® M a |A), 



(3.3) 



arising from the splitting of the algebra generated by the S Q 's: A — Q) Je jI nj ® M.(dj,C). This splitting of the 
algebra has allowed both DFSs and QECCs to be put into a similar framework J39||l2|]. We will now show that the 
DFS condition on the algebra A generated by the S Q can be converted into a stabilizer condition on the complex Lie 
algebra generated by the S Q 's. We define the continuous DFS-stabilizcr D(z7) as 



■V N 



exp 



E 



v a (S a - 1 <g> M a ) 



v a G C. 



Clearly, if the DFS condition Eq. ( |3~5| ) is fulfilled for a set of states \fi) ® |A), then 

D(w)|/x) ® |A) = ® | A) Wv a G v. 



(3.4) 



(3.5) 



Thus the DFS condition implies that the D(t7) stabilize the DFS. Further if E q. (3.5) holds then in particular it must 
hold for a v which has only one non-vanishing component vp. Thus Eq. (3J5) implies that D(0, . . . , 0, vp, . . . 0)|/z) €5 
| A) = |/z)(S>|A). Recalling that exp[-] is a one-to-one mapping from a neighborhood of the zer o matrix to a neighborhood 
of the identity matrix, it follows that there must exist a small enough v a such that Eq. (3.5) implies the DFS condition 
Eq. (3.3). Thus we see that we can convert the DFS condition into a condition on the stabilizer of the complex Lie 
algebra generated by the S a — I ® M Q 's: 



|*)GDFS iff D(v)|tf) = I*) VwGC 



n 



(3.6) 



In some cases we will be able to pick a finite subgroup from elements of T)(v) which constitutes a stabilizer. We 
will mention these instances in the following sections. However, apart from the conceptual framework, our main 
motivation to introduce the stabilizer for a DFS is to be able to analyze the errors which a DFS (i) detects /cor rect s 
(as a QECC), and (ii) those which it avoids (passive error correction). The continuous stabilizer provided in Eq. (3.4) 
will be sufficient to study these errors. 

As in the previous paragraph, errors E Q (i) that anticommutc with an clement in the stabilizer will take codewords 
to subspaces that are orthogonal to the code. These errors will be detectable (and correctable if eJjE q anticommutes 
with a stabilizer element) ||. In order to identify the QECC-properties of a DFS, it will be convenient to look for 
elements of the Pauli group among the D(w). 

A code C with stabilizer D(w) will avoid errors of type (ii) in its stabilizer in the sense that, if all of the Kraus 
operators of a given decoherence process can be expanded over stabilizer elements Aj(f) = J^,n ei{v, t)D(v)dv, then 



p(t) = ^A i (i)p(0)A](t)=^ / 



(v, t)dv 



p(0). 



(3.7) 



8 



The normalization condition ^ i Aj(t)t Aj(i) = I then implies that ^ 



= 1. Consequently, as expected, 



/ ei(v,t)dv 

the DFS does not evolve. Hence we see that the stabilizer provides an efficient method for identifying the errors which 



a code avoids. In later sections we analyze the concrete form of the stabilizer Eq. (3.4) for the error models studied 
in this paper. 

IV. THE COMMUTANT AND UNIVERSAL QUANTUM COMPUTATION ON A DFS 

A DFS is a promising way to store quantum information in a robust fashion pij . From the perspective of quantum 
computation however, it is even more important to be able to controllably transform states in a DFS, if it is to be 
truly useful for quantum information processing. More specifically, to perform quantum algorithms on a DFS one 
has to be able to perform universal quantum computation using decoherence-free states. The notion of universal 
computation is the following: with a restricted set of operations or interactions at hand, one wishes to implement any 
unitary transformation on the given Hilbert space, to an arbitrary degree of accuracy. From a physical implementation 
perspective it seems clear that the operations used (gates) should be limited to at most two-body interactions. In 
particular we wish to identify a finite set of such gates that is universal on a DFS. 

Since we do not wish to implement active QECC, we impose a very stringent requirement on the operations we 
shall allow for computation using DFSs: we do not allow gates that ever take the decoherence free states outside the 
DFS, where the states would decohere under the noise- process considered.^ As a first step towards this goal we thus 
need to be able to identify the physical operations which perform transformations entirely within the DFS. 

A. Operations that Preserve the DFS 

There are essentially two equivalent approaches to identify the "encoded" operations that preserve a DFS. One is via 
the normalizer of the stabilizer of a code |p^| ; the second is via the commutant of the f-closed algebra generated by 
the error operators ]39| . Both will be briefly reviewed here. 

Computation on a stabilizer DFS: The stabilizer formalism is very useful for identifying allowed gates that take 
codewords to codewords p3|. An operation U keeps code- words inside the code-space, if and only if the trans- 



formed state U|\I/) is an element of the code C. Thus, using the stabilizer condition (3.1) for codes with stabilizer S 
and C = : S|tf) = |W) VS E S}, we have 

U|*)eC iff SU|#)=U|#) VSeS. (4.1) 

This implies U _1 SU|5') = \9) and so U _1 SU e S: Allowed operations U transform stabilizer elements S by 
conjugation into stabilizer elements; U is in the normalizer of S (if S is a group). If we restrict the allowed operations 
to gates in the Pauli-group (as is done in 0|), then the allowed gates U will fix the stabilizer pointwise (element by 
element). In the case of DFS with a continuous stabilizer D(u), the above translates to the following condition p2| 

UD^U^D^p)), (4.2) 

together with the requirement that T)(v' '(v)) must cover S. To satisfy the covering condition, it is sufficient to have 
^{v) be a one-to-one mapping. 

Eq. ( |4.2j ), derived by generalizing concepts from the theory of stabilizers in the Pauli group, is a condition that 
allows one to identify gates U that transform codewords to codewords. In a physical implementation these gates will 
be realized by turning on Hamiltonians H between physical qubits for a certain time t: U(t) = e" JtH . So far we 
only required that the action of the gate preserve the subspace at the conclusion of the gate operation, but not that 
the subspace be preserved throughout the entire duration of the gate operation. The stabilizer approach allows us to 
further identify the more restrictive set of Hamiltonians that keep the states within the DFS throughout the entire 
switching time of the gate. As a result, in the limit of ideal gates, the entire system is free from noise at all times. 
This is different from QECC, since there errors continuously take the codewords outside of the code-space Q, and 
hence error correction needs to be applied frequently even in the limit of perfect gate operations. Imperfections in 



2 We shall lift this requirement in Sec. VIII 
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gate operations can be dealt with in the DFS approach by concatenation with a QECC p2|, as shown explicitly for 
the exchange interaction injS]. 



By rewriting condition (4.2) as U(<)D(u) = T)(v'(v, t))U(t), taking the derivative with respect to t and evaluating 
at t = wc obtain HD(w) = T)(v'(v, 0))H + i ( so that: 

Theorem 1 — A sufficient condition for the generating Hamiltonian to keep a state at all times entirely within the 
DFS is HD(?) = D(i? '(w))H where v'{v) is one-to-one and time- independent. 

For most applications we will only need gates that commute with all stabilizer elements. The condition for the 
generating Hamiltonian then simplifies to HD(v) = D(w)H. 

Computation on irreducible subspaces: We can derive conditions to identify allowed gates on a DFS by using the 
representation theoretic approach developed in |59) , and section H. Recall that the decomposition of the algebra 
A = 0j e j l nj ®M.{dj, C) generated by the errors {So,} induces a splitting of the Hilbert space Hs = J2jej C" J ®C J 
into subspaces possessing a tensor product structure suitable to isolate decoherence-free subsystems. The set of 
operators in the commutant of A, A' = {X : [X, A] = 0,VA 6 A} = A' = (§)j e jM(nj,€) ® ldj, obviously generate 
transformations that affect the codespace only. In particular, they take states in a DFS to other states in that same 
DFS. A' is generated by operators which commute with the S Q . Again, our goal is to find gates that act within a 
DFS during the entire switching time, and to this end we need to identify Hermitian operators H in A' to generate 
an evolution U(t) = exp[— itH] on the DFS. 

Theorem 2 — A sufficient condition for a Hamiltonian H to generate dynamics U(t) = exp[— itH] which preserves 
a DFS is that H be in the commutant of the algebra A. 

However, because we can only use one particular DFS (corresponding to a specific K € Sf) to store quantum 
information (the coherences between superpositions of different DFSs are not protected) , the operators which commute 
with the S Q 's are not the only operators which perform non-trivial operations on a specific DFS. The operations in A' 
preserve all DFSs in parallel. However, if we restrict our system to only one such DFS, we do not need any constraints 
on the evolution of the other subspaces. It is then possible to construct a necessary and sufficient condition for a 
Hamiltonian by modifying the commutant to: 

T = (M(n K) €)®I dK )®M(d-d K n K ,€) (4.3) 

where dx = dim(7is) and just leaves the specific DFS (K) invariant. 

Theorem 3 — A necessary and sufficient condition for a Hamiltonian H to generate dynamics which preserves a 
DFS corresponding to the irreducible representation K S J , is H G T. 

We will use both the stabilizer and the commutant approaches, to find a set of universal gates for decoherence 
processes of physical relevance. In the cases discussed in this paper, any one of the two approaches is clearly sufficient 
and we do not need all theorems in full generality. However we provide here a general framework and the tools 
required to analyze DFS and QECC stemming from any error model. 

Finally we should point out again that from a practical perspective, it is crucial to look for the Hermitian operations 
which perform nontrivial operations on the DFS and which correspond to only one or two-body physical interactions. 



Without this requirement, it is clear that one can always 29 construct a set of Hamiltonians (satisfying the conditions 



of Theorem 2) which span the allowed operations on a DFS. A primary goal of this paper is therefore to construct such 
one and two-body Hamiltonians for specific decoherence mechanisms, in order to achieve true universal computation 
on the corresponding DFSs. 



B. Universality and Composition of Allowed Operations 

Using the tools developed in the previous subsection, we can now find local one-and-two qubit gates that represent 
encoded operations on DFSs. However, in general, a discrete set of gates applied in alternation is not sufficient to 
generate a universal set of gates. Nor is it sufficient to obtain every encoded unitary operation exactly. Furthermore, 
for analysis of the complexity of computations performed with a given universal set of gates, it is essential to keep 
under control the number of operations needed to achieve a certain gate within a desired accuracy. In the theory of 
universality (e.g., ||(J) the composition laws of operations have been analyzed extensively. We will review the essential 
results relevant for our purposes here. 

Let us assume that we have a set of (up to two-body) Hamiltonians H = {H^ : i = 1, . . . , M} that take DFS states 
to DFS states. We will construct gates using the following composition laws: 

1. Arbitrary phases: Any interaction can be switched on for an arbitrary time. Thus any gate of the form 
U(i) = cxp (— iiH-i) can be implemented. 
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2. Trotter formula: Gates performing sums of Hamiltonians are implemented by using the short-time approxi- 
mation to the Trotter formula exp [i(tiHi + t2Hj)]= lim„_ s . 0O [exp (i&HA exp u-^Hjj] : 

e i(tiH s +t2H 3 )/n _ e itiHi/n e it 2 H 3 /n _|_ Qf_\ U 

V? 

This is achieved by quickly turning on and off the two interactions FT, with appropriate ratios of duration 
times. An alternative, direct, way of implementing this gate is to switch on the two interactions simultaneously 
for the appropriate time-intervals. 

3. Commutator: It is possible to implement the commutator of operations that are already achievable. This is 
a consequence of the Lie product formula 

cxp[Hi,Hj] = lim [exp (itli/y/n) exp (iHj/y/nj exp (— iHi/y/n) exp (— iHj/v^)]" , 
which has the short-time approximation 

gtlHi^Hj-J/n^gitHi/v^gitHj/v^g-itHi/v^Tg-itHj/v^r _|_ qi \ \ /^gN 

Again, the gate e lt ^~^ Hi ' Hj ^ can be implemented to high precision by alternately switching on and off the 
appropriate two interactions with a specific duration ratio.^ 

4. Conjugation by unitary evolution: Another useful action in constructing universal sets of gates comes from 
the observation that if a specific gate U and its inverse ut can be implemented, then any Hamiltonian H which 
can be implemented can be modified by performing U before, and after the gate exp(— itH). This gives rise 
to the transformed Hamiltonian 

Uexp(-i<H)U t = exp(-i<UHU t ) = cxp(-itH off ). (4.6) 



Note that the laws (1-3) correspond to closing the set of allowed Hamiltonians as a Lie-algebra (scalar multiplication, 
addition and Lie-commutators can be obtained out of the given Hamiltonians) . 

If (a subset of) the composition laws (1-4) acting on the set H give rise to a set of gates that is dense in the group 
SU(cIk) (via successive application of these gates), where oIk is the dimension of the DFS, then we shall refer to H as 
a universal set of generators. Equivalently, this means that H generates the Lie-algebra su(d,K ) (traceless matrices) 
via scalar multiplication, addition, Lie- commutators, and conjugation by unitaries. The generators of this algebra can 
be obtained from H by these operations. 

For all practical applications and implementations of algorithms, we will only be interested to approximate a certain 
gate sequence with a given accuracy. Note that the composition laws (2) and (3) use only repeated applications of (1) 
in order to approximate a certain gate. We can replace the requirement to perform an arbitrary phase, (1), by noting 
that e** Hi is dense in the group {exp(iiHi) : t S [0, 2tt)} if t is an irrational multiple of tt. Repeated application of 
that gate can then approximate an arbitrary phase to any desired accuracy. Thus we can in principle restrict our 
available gates to {exp(itjHj)}, together with fixed irrational switching times U. Repeated application of these gates 
can then be used to approximate any operation in SU(dK) to arbitrary accuracy. 

In order to prove that a set H generates a universal set of Hamiltonians, we use the fact that a large group of 
universal sets have already been identified [ p5| , p6|j5l| ] . It suffices to show that H generates one of these sets, in order 
to prove that H is a universal set of generators. We will use the fact that the set of one qubit operations SU(2) is 
generated by any two arbitrary rotations with irrational phase, around two non-parallel axes. Alternatively, if we are 
given these two rotations with any phase, then an Euler-angle construction can be used to yield any gate in SU (2) 
by application of a small number of rotations (three if the axes are orthogonal). In addition we will use (and prove) a 
lemma {Enlarging Lemma, Appendix ^) that allows extension to su(n + l) of a given su(n) acting on an n-dimensional 
subspace of a Hilbert space of dimension n + 1, with the help of an additional su(2). 



3 Note that in order to implement e * tA we would use e li>A = 7 and implement e*^ *' A instead. This depends on A having 
rationally related eigenvalues, which will always be the case for the Hamiltonians of interest to us. 
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In order to use this approach to universality, it is crucial to have bounds on the length of the gate sequences 
approximating a certain gate in terms of the desired accuracy. This is all the more important if one universal set is to 
be replaced by any other with only polynomial overhead in the number of gates applied, for otherwise the complexity 
classes would not be robust under the exchange of one set for another. The whole notion of universality would 
then by questionable. The following key theorem proved independently by Solovay and Kitaev (see |M) establishes 
the equivalence of universal sets, and provides bounds on the length of gate sequences for a desired accuracy of 
approximation. In order to quantify the accuracy of an approximation, we need to define a distance on matrices. 
Since our matrices act in a space of given (finite) dimension dx , any metric is as good as any other. For example, we 

can use the trace-norm d(U, V) = y^l — ^rRe [Tr(UtV)] . A matrix V is then said to approximate a transformation 

U to accuracy e if rf(U, V) < e. 

Theorem (Solovay-Kitaev) — Given a set of gates that is dense in SU(2 k ) and closed under Hermitian conjugation, 
any gate U in SU(2 k ) can be approximated to an accuracy e with a sequence of poly [log(l/e)] gates from the set. 

DFS- Corollary to the Solovay-Kitaev Theorem — Assume that the DFS encodes a d^-dimensional system into n 
physical qubits. Given that one can exactly implement the gate set {e ltiHi : Hj S H}, [£< are (fixed) irrational 
multiples of n, and H is a universal generating set] it is possible to approximate any gate in SU(di() (any encoded 
operation) using m = poly [log(l/e)] gates. 

Furthermore, if we can only implement the given gates approximately, say to an accuracy 5, we will still be able 
to approximate the target gate: It is known that a sequence of m J-imprecise unitary matrices is (in some norm) at 
most distance mS far from the desired gate. If a sequence of exactly implemented gates Ui , . . . , U m approximates 
a target gate U up to e, and instead of Ui, . . . , U m , we use gates that are at most some distance S apart, then the 
total sequence will be at most e + m8 = e + poly [log(l/e)] <5 apart from U. If we make sure that S < epoly [log(l/e)] 
then the <5-faulty sequence will still approximate U to a precision 2e. 

If we further assume that the physical interaction that we switch on and off is given by the device and is unlikely 
to change its form, then the imprecision of the gate comes entirely through the coupling strength and the interaction 
time, i.e. a faulty gate is of the form TJf = e I (0+ A 0) H ) where U = e^ 11 is the unperturbed gate. The distance 

d(U,U/) = J 1 - Re [Tr(e iA ^ H )] = J 1 - -j- Re [Tr(cos A01 + i sin A0H)] 
V d K V "if 

= y/l -cosA0 = V2sin(A< ? !)/2) « ^ = 5, (4.7) 

v2 

is proportional to the error A</> of the product of coupling-strength and interaction time. This translates to (nearly) 
linear behavior in the desired final accuracy e. 



V. COLLECTIVE DECOHERENCE 



We now focus on a particularly interesting and useful model of a DFS. This is the case of collective decoherence 
on n qubits. We distinguish between two forms of collective decoherence. The first, and simpler, type of collective 
decoherence is weak collective decoherence (WCD). We define the collective operators as 

n 

j'=i 

where o~ J a denotes a tensor product of the a th Pauli matrix, a — x,y, z, 

a * = (°il) ff "=(! o ! ) = ( -1 ) ( 5 - 2 ) 

(in the basis spanned by a z eigenstates |0) and |1)) operating on the j th qubit, and the identity on all of the other 
qubits. WCD is the situation in which only one collective operator S a is involved in the coupling to the bath, i.e., 
H/ = S a ® B. 

The second, more general type of collective decoherence is strong collective decoherence (SCD). We define SCD 
as the general situation in which the interaction Hamiltonian is given by H/ = X) Q =i ® B a . The S a provide a 
representation of the Lie algebra su(2): 

[S Q , Sp] — — 2ie Q/ 3 7 S 7 (5.3) 
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The B Q 's are not required to be linearly independent. 

Both of these collective decoherence mechanisms are expected to arise from the physical condition that the bath 



cannot distinguish the system qubits 13 14JlaJ21]. If there are n qubits interacting with a bath, the most general 



interaction Hamiltonian linear in the a l a is given by 

n 

H/=£ E <® B i,<* ( 5 - 4 ) 

i— 1 a—x,y,z 

where the Bi jQ are bath operators. If the bath cannot distinguish between the system qubits, then Bi jQ should not 
depend on i and the Hamiltonian becomes H/ = J2 a =x y z ^ a ® ^o> i- e -> strong collective decoherence. 

As a concrete example of such collective decoherence, consider the situation in which the bath is the electromagnetic 
field, and the wavelength of the transition between the states of the qubits is larger than the spacing between the 
qubits. The electromagnetic field will interact with each of these qubits in an identical manner, because the field 
strength over a single wavelength will not vary substantially. This gives rise to the well- known phenomena of Dickc 
super- and sub-radiance |5^| . Whenever the bath is a field whose energy is dependent on its wavelength and this 
wavelength is much greater than the spacing between the qubits, one should expect collective decoherence to be 
the dominant decoherence mechanism. It is natural to expect this to be the case for condensed-phase high-purity 
materials at low temperatures. However, to the best of our knowledge at present a rigorous study quantifying the 



relevant parameter ranges for this interesting condition to hold in specific materials is still lacking (see Refs. |E7 28 
for an application to quantum dots, though). 

VI. THE ABELIAN CASE: WEAK COLLECTIVE DECOHERENCE 

For a decoherence mechanism with only one operator S a coupling to the bath, the implementation and discussion of 
universal computation with local interactions is simpler than in the general case, because we can work in the basis that 
diagonalizes S Q (S a is necessarily Hermitian in the Hamiltonian model we consider here). The algebra generated by 
S Q is abelian and reduces to one-dimensional (irreducible) subalgebras corresponding to the eigenvalues of S a . More 
specifically, Ai = ® A/ l nj <8 M{Xj), where Xj is the J th eigenvalue with degeneracy nj, and M(\j) is the algebra 
generated by Xj. A4(Xj) acts by multiplying the corresponding vector by Xj. In this situation the DF subsystems 
are only of the DF subspace type. This simpler case of weak collective decoherence allows us to present a treatment 
with examples, that will make the general case of strong collective decoherence (SCD) more intuitive. 

In the following we will, without loss of generality, focus on the case S a = S 2 = X)l-=i °z -0 This operator is already 
diagonal in the computational basis (the eigenstates are bitstrings of qubits in either |0) or |1)). Since a k z acting on 
the k th qubit contributes 1 if the qubit is |0), and —1 if the qubit is |1), the eigenvalue of a bitstring is #0 — #1 (the 
number of zeroes minus the number of ones), and the eigenvalues of S z are {n, n — 2, . . . , — n + 2, — n}. 

The degeneracy nj of the eigenspace corresponding to an eigenvalue 

Xj = n-2J (6.1) 



•»-{') < 6 - 2 > 

(the number of different bitstrings with n — J zeroes and J ones) . The abelian algebra generated by S z thus splits into 
one-dimensional subalgebras with degeneracy nj. The largest decoherence- free subspaces in this situation correspond 
to the space spanned by bitstring-vectors where the number of zeroes and the number of ones are either the same (n 
even), or differ by one (n odd). 

A. The Stabilizer and Error Correction Properties 



Following the formalism developed in Section (III) we find, using Eq. (3.4) with v = i6 (9 can be complex), the 



stabilizer for the weak case corresponding to a DFS with eigenvalue Xj to be 



The cases a = x (y) follow by applying a bitwise Hadamard (Hadamard+phase) transform to the code. 
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Z® n (0) = exp[^(S, - Xjl)} =(g)e-^ je (Icos6 + <7 k z ishi9) = e - iXje V{6)® n (6.3) 

fc=i 

where 

nO) = ( ^ e i ) . (6.4) 

For strictly real 9 some of the errors which are protected against are simply collective rotations about the a z axis 
(and an irrelevant global phase). For strictly imaginary 9 we find that the errors which are protected against are 
contracting collective errors of the form diag(e e , e~ e ), i.e., they result in loss of norm of the wave function. Any 
physical process with Kraus operators that are linear combinations of these errors will therefore not affect the DFS. 

This is the right framework in which to present another form of the stabilizer. We note that in the case of weak 
collective decoherence, we can find a stabilizer group with a finite number of elements. Define 

„ (2m \ /exp(^f) \ ,„ 

Z i= e ^{ — ^)-{ exp(-^)J- ^ 

Then the n-element group Z n generated by exp (— i2n\j/n) Zf^ is a stabilizer for the DFS corresponding to the 
eigenvalue A,/. To see that 

^{-*^\zt n \*) = \V) iff |*) € DFS(A J ), (6.6) 

note that a Z\i n acting on a |0) contributes exp(27ri/n) to the total phase, whereas Z\j n acting on a |1) contributes 
exp(— 2m/n), So Zfj l n gives a total phase of exp (27n(#|0) — #|l))/n) = exp(27riAj/n) when acting on a bitstring. 
This stabilizer and Eq. ( |6.6| ) provide a simple criterion to check whether a state is in a DFS or not. 

Let us now briefly comment on the error-correction and detection properties of the code in the WCD case. The 
stabilizer elements are all diagonal, and equal to a tensor product of identical 1-qubit operators. The element Z® n is 
in the stabilizer and anticommutes with odd-number X and Y errors. So odd-number qubit bit-flips are detectable 
errors. However the code is not able to detect any form of error involving Z's and even-number X's and Y's, since 
any such error commutes with all elements in the stabilizer. 

B. Nontrivial Operations 

Observe that the algebra A in the WCD case is generated entirely by S z . Hence, by Theorem 2, the DFS-preserving 
operations are those that are in the commutant of S z . For single-body Hamiltonians it is easy to see that the only 
non-trivial such set is formed by interactions proportional to o\ operators. As for two-qubit Hamiltonians, it is simpler 



to use Theorem 1 and the expression (6.3) for the stabilizer. We are then looking for 4x4 Hermitian matrices that 
commute with P(#)® 2 ; these are of the form 



Ty(zi, z 2 , z 3 , z 4 , h) 



(z x 











\ 





Z2 


h 










h* 


Z3 







V o 










) 



(6.7) 



where Ty acts on qubits i and j only. Here Zi is real, h is complex, and the row space is spanned by the i th and j th 
qubit basis {|00), |01), |10), |11)}. We note that systems with an internal Hamiltonian of the Hciscnbcrg type, 



n 



1 



H Hc is = > tjo-l + - > v Jyai ■ <Jj , (6.8) 



/ j J~ Z • / j 

3=1 " i,3=l 



have exactly the correct form for any pair of spins Indeed, it is not hard to see that [HHeis,S 2 ] = |25|. The 
Heisenberg Hamiltonian is ubiquitous, and appears, e.g., in NMR. This means that the natural evolution of NMR 
systems under WCD preserves the DFS, and implements a non-trivial computation. 
The specific case 
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E; 



T« (1,0, 0,1,1) 



( 1 
10 
10 

Vo o o i 



(6.9) 



which flips the two states 1 01) and 1 10} of qubits i and j and leaves the other two states invariant, is especially 
important: it is the exchange interaction. The other interactions we employ are 



T£ =Ty(l,0,0,0,0) = diag(l, 0,0,0) 



Tg = T«(0,0,0,l,0) = diag(0, 0,0,1), 
which introduce a phase on the state 1 00) (P) and 1 1 1} (Q) of qubits i and j; and 

Zia = T 12 (0, 0, 1,0, 0)} = diag(0, 0, 1, 0). 



(6.10) 



(6.11) 



In the following we show that these special interactions are sufficient to obtain a universal generating set operating 
entirely within a weak-collective DFS. 



C. Universal Quantum Computation inside the Weak-Collective DFS 

Let DFS n (K) denote the DFS on n physical qubits with eigenvalue K. We show here that 

H = {E i , i+ i,T^ fl ,T« +1 :i = l,...,n-l,Zi 2 } (6.12) 

is a universal generating set for any of the DFSs occurring in a system of n physical qubits. It is convenient to 
work directly with the Hamiltonians, and to show that H gives rise to the Lie-algebra su(cIk) on each DFS„(injvia 
scalar multiplication, addition, and Lie-commutator; see the allowed compositions of operations (1-3) in Section |VB 



Exponentiation then gives the group SU(oIk) on the DFS. We will proceed by induction on n, the number of physical 
qubits, building the DFS-states of n qubits out of DFS-states for n — 1 qubits. A graphical rep resen tation of this 
construction is useful (and will also generalize to the strong case presented in the following section VII): see Fig. (Q). 



We have seen that in the WCD case the DFS states are simply bitstrings of n qubits in either |0) or |1). The 
different n-qubit DFSs are labeled by their eigenvalue 

A, 7 = #0-#1 = J n . (6.13) 

To obtain a DFS-state of n qubits out of a DFS-state of n — 1 qubits corresponding to J n _i we can either add the 
n th qubit as |0) (J n = J n -i + 1) or as |1) (J n = J n -i — 1). Each DFS-state can be built sequentially from the first 
qubit onward by adding successively JO) or |1), and is uniquely defined by a sequence Ji, .. ., J n of eigenvalues. In 
the graphical representation of Fig. (pi) the horizontal axis marks n, the number of qubits up to which the state is 
already built, and the vertical axis shows J n , the difference #0 — #1 up to the n th qubit. Adding a |0) at the n + 1 
step will correspond to a line pointing upwards, adding a |1) to a line pointing down. Each DFS-state of nqubits with 
eigenvalue Xj — J n is thus in one-to-one correspondence with a path on the lattice from the origin to in, J n ). 

Consider the first non-trivial case, n = 2, which gives rise to one DFS-qubit: DFS2(0). This corresponds to the two 
states |0l) = 1 01) [path 2 in Fig. (|l|)] and = 1 10) (path 3) with J 2 = 0. The remaining Hilbert space is spanned 
by the one-dimensional DFS2(2) |00) (path 1) corresponding to J2 = 2, and DFS2(— 2) 1 1 1) (path 4) corresponding 
to J2 = —2. The exchange E12 flips \0l) and (path 2 and 3), and leaves the other two paths unchanged. The 
interaction A 12 = diag(0,0, 1,0) induces a phase on |1^) = 1 10) (path 3). Their commutator forms an encoded a y 
acting entirely within the DFS2(0) subspace. Its commutator with E12 in turn forms an encoded a z with the same 
property. Together they form the (encoded) Lie algebra su(2) acting entirely within this DFS. The Lie algebra is 
completed by forming the commutator between these Y and Z operations. To summarize: 




Y 12 = i[A,E ia ]= " " * " (6.14) 



Z12 = i[Ei2, Y12] 

Xi 2 = i[Y 12 ,Z 12 ] (6.15) 
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We call the property of acting entirely within the specified DFS independence, meaning that the corresponding Hamil- 
tonian has zero entries in the rows and columns corresponding to the other DFSs [DFS2(2)=|00) and DFS2(— 2)=|11) 
in this case]. When the Hamiltonian is exponentiated, the corresponding gate will act as identity on all DFSs except 
DFS 2 (0). 

To summarize these considerations, the Lie-algebra formed by Hg = {X, Z} is su(2), and generates SU(2) on 
DFS2(0) by exponentiation. In addition, this is an independent SU(2), namely, these operations act as identity on 
the other DFSs: when written as matrices over the basis of DFS-states, their generators in Hq have zeroes in the rows 
and columns corresponding to all other DFSs. 

In the following we show how this construction generalizes to n > 2 qubits, by proving the following theorem: 

Theorem 4 — For any n > 2 qubits undergoing weak collective decoherence, there exist sets of Hamiltonians H™ n 
[obtained from H of Eq. (6.12) via scalar multiplication, addition, and Lie-commutator] acting as su(dj n ) on the DFS 
corresponding to the eigenvalue J n . Furthermore each set acts independently on this DFS only (i.e., with zeroes in 
the matrix representation corresponding to their action on the other DFSs). 

Before proving this theorem, we first explain in detail the steps taken in order to go from the n = 2 to the n = 3 
case, so as to make the general induction procedure more transparent. 

The structure of the DFSs for n = 2 and 3 qubits is: 



DFS 2 (2) = {|00)}, DFS 2 (0) = 



DFS 3 (3) = {|000)}, DFS 3 (1) 



|01) 
1 10) 



DFS 2 (-2) = {|11)} 



looi) r ion) 

|010) , DFS 3 (-1) = { |101) , DFS 3 (-3) = {|111)}. 
|100) I |110) 



(6.16) 



DFS 3 (3) is obtained by appending a |0) to DFS 2 (2). Similarly DFS 3 (-3) is obtained by appending a |1) to DFS 2 (-2). 
Graphically, this corresponds to moving along the only allowed pathway from DFS2(2) [DFS 2 (— 2)] to DFS 3 (3) 
[DFS 3 (— 3)], as shown in Fig. (jl]). The lowest and highest Aj for n qubits will always be made up of the single 
pathway connecting the lowest and highest Aj for n — 1 qubits. The structure of DFS 3 (±1) is only slightly more 
complicated. DFS 3 (1) is made up of one state, 1 001) , which comes from appending a |1) (moving down) to DFS 2 (2). 
We call 1 001) a "Top-state" in DFS 3 (1). The two other states, 1 010) and 1 100) , come from appending |0) (moving up) 
to DFS 2 (0). Similarly, we call |010) and 1 100) "Bottom-states" in DFS 3 (1). DFS 3 (— 1) is constructed in an analogous 
manner (Fig. |l|). 

We showed above that it is possible to perform independent su(2) operations on DFS 2 (0). DFS2(±2) are also 
both acted upon independently, but because they are one-dimensional subspaces, independence implies that su(2) 
operations annihilate them. Since the states {|010>, |100)> € DFS 3 (1) and the states {|011), |101)} € DFS 3 (-1) 
both have {|01), |10)} £ DFS 2 (0) as their first two qubits, one immediate consequence of the independent action on 
DFS 2 (0) is that one can simultaneously perform su(2) operations on the corresponding daughter subspaces created 
by expanding DFS 2 (0) into DFS 3 (±1). The first step in the general inductive proof is to eliminate this simultaneous 
action, and to act independently on each of these subspaces (the "independence step"). To see how this is achieved, 
it is convenient to represent the operators acting on the 8-dimensional Hilbert space of 3 qubits in the basis of the 4 
DFSs: 



000 


001 010 100|011 101 110 


111 


M 3 










Mi 










M-i 










M_ 3 



The simultaneous action on DFS 3 (±1) can now be visualized in terms of both M±\ being non-zero. Let us show 
how to obtain an action where, say, just Mi is non-zero. This can be achieved by applying the commutator of two 
operators with the property that their intersection has non- vanishing action just on M\. This is true for the 
and X12 Hamiltonians: annihilates every state except those that are |00) over qubits 2 and 3, namely 1 100) € 
DFS 3 (1) and 1 000) €DFS 3 (3). This implies that the only non-zero blocks in its matrix are 



M 3 (Tf 3 ) = l, M 1 (Tg) 




(6.17) 
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On the other hand, Xi 2 is non-zero only on those states that are 1 01} or 1 10) on qubits 1 and 2. Therefore it will 
be non-zero on all 3-qubit states that have 1 01} or 1 10) as "parents". This means that in its matrix representation 
M ± 3 = and 



Mi(Xi 2 ) = 1, M_!(X 12 ) =10 . (6.18) 





Clearly, taking the product of T^ 3 and X12 leaves non-zero just the lower 2x2 block of Mi, and this is the crucial 
point: it shows that an independent action on DFS3(1) can be obtained by forming their commutator. Specifically, 
since the lower 2x2 block of Mi(T% 3 ) is just § (I - a z ): 

*[ T 23' X 12] = Y{| 10 0),|010>}> (6.19) 

i.e., this commutator acts as an encoded a y inside the { 1 100) , |010)} subspace of DFS3(1). Similarly, Z{| 100 ) i |oio)} = 
|[Y{|ioo},|oio»,Xi 2 ]. Together {Y { | 100)i | 010)} , Z { , 10 o> , |oio> } } generate su(2) acting independently on the {|100), |010)} 
subspace of DFS3(1), which we achieved by subtracting out the action on DFSa(— 1). 

In an analogous manner, an independent su(2) can be produced on the {|011), |101)} subspace of DFS3(— 1) by 
using the Hamiltonians acting on DFS2(0) in conjunction with T^ 3 to subtract out the su(2) action on DFS3(1).^| 
Thus we can obtain independent action for each of the daughters of DFS2(0), i.e., separate actions on the subspace 
spanned by {|010), |100)} and {|011), |101)}. 

Having established independent action on the two subspaces of DFSs(l) and DFSs(— 1) arising from DFS2(0), we 
need only show that we can obtain the full action on DFSa(l) and DFS3(— 1). For DFS3(1) we need to mix the 
subspace {|010), 1 100) } over which we can already perform independent su(2), with the 1 001) state. To do so, note 
that the effect of the exchange operation E23 is to flip |001) and 1 010) , and leave 1 100) invariant. I.e., the matrix 
representation of E23 is 



Mx(E 23 ) = 1 . (6.20) 




Unfortunately, E23 has a simultaneous action on DFSa(— 1). This, however, is not a problem, since we have already 
constructed an independent su(2) on DFSa(l) elements. Thus we can eliminate the simultaneous action by simply 
forming commutators with these su(2) elements. The Lie algebra generated by these commutators will act indepen- 
dently on all of DFSa(l). In fact we claim this Lie algebra to be all of su(3) (see Appendix for a general proof). 
In other words, the Lie algebra spanned by the su(2) elements {a x , a y , a z } acting on the subspace { 1 100) , |010)}, to- 
gether with the exchange operation E23, generate all of su(3) independently on DFS(l). A similar argument holds for 
DFSa(— 1). This construction illustrates the induction step: we have shown that it is possible to perform independent 
su(dK) actions on all four of the DFSs(A) (K — ±3, ±1), given that we can perform independent action on the three 
DFS 2 (A) (K = ±1,0). In Fig. (|J) we have further illustrated these considerations by depicting the action of exchange 
on two the 4-qubit DFSs. Let us now proceed to the general proof. 
Proof — By induction. 

The case n — 2 already treated above will serve to initialize the induction. Assume now that the theorem is true 
for n — 1 qubits and let us show that it is then true for n qubits as well. 

First note that each DFS„(A) is constructed either from the DFS„_i(A" — 1) (to its lower left) by adding a |0) for 
the n th qubit, or from DFS„_i(A + 1) (to its upper left) by adding a |1): the states in DFS„(A) correspond to all 
paths ending in (n, K) that either come from below (B) or from the top (T). See Fig. (0). 



5 Since T23 annihilates every state except those that are |11) over qubits 2 and 3, namely |011) £ DFSa( — 1) and 
1 1 1 1) GDFSs(— 3), the only non-zero blocks in its matrix are 

M_ 3 (Tg,) = l, M_i(T£) = 




17 



If we apply a certain gate U = exp(iHi) to DFS„_i (K+l), then this operation will induce the same U on DFS„(A), 
by acting on all paths (states) entering DFS„(A) from above. At the same time U is induced on DFS„(A + 2) by 
acting on all paths entering this DFS from below. So, U affects two DFSs simultaneously. In other words, the set of 
valid Hamiltonians H^~7 \ [acting onra-1 qubits and generating su(cIk+i)] on DFS„_i(A + 1), that we are given by 
the induction hypothesis, induces a simultaneous action of su(cIk+i) on DFS„(A) (on the paths coming from above 
only) and DFS n (A + 2) (on the paths coming from below only). Additionally, it does not affect any other n-qubit 
DFS, since we assumed that the action on DFS n -i(A + 1) was independent, and the only n-qubit DFSs built from 
DFS n _i(A + 1) are DFS„(A) and DFS„(A + 2). These considerations are depicted schematically in Fig. (||). 

We now show how to annihilate, for a given non-trivial (i.e., dimension > 1) DFS„(A), the unwanted simultaneous 
action on other DFSs (the "independence step"). Then we proceed to obtain the full sit(dft-), by using the su(dK±i) 
on DFS„_i(A± 1) that are given by the induction hypothesis (the "mixing step"). 

1. Independence 

Let us call all the tx paths converging on DFS„(A) from above "Top-states", or T-states for short, and the bj< 
paths converging from below "Bottom- (or B) states" (recall that there is a 1-to-l correspondence between paths and 
states). The total number of paths converging on a given DFS is exactly its dimension, so c?a" = tx + ok- By using 
the induction hypothesis on DFS„_i(A" + 1) we can obtain su{tn) (gener ated by H™^\) on the T-states of DFS n (A), 
which will simultaneously affect the B-states in the higher lying DFS„(A + 2) as su(ok+2) (note that tjc = &a"+2)- 
The set H^rA is non-empty only if n — 3 > K + 1 > — (n — 3) [because the "highest" and "lowest" DFS are always 
one-dimensional and su(l) — 0). If this holds then DFS„(A + 2) "above" DFS„(A) is non-trivial (dimension > 1), 
and there are paths in DFS„(A) ending in |11) ("down, down"). This is exactly the situation in which we can use 
Tn-i n t° wipe out the unwanted action on DFS„(A + 2): recall that T®_ 1 n annihilates all states except those ending 
in |11), and therefore affects non-trivially only these special T-states in each DFS. Since the operations in H^r^ affect 
only B-states on DFS„(A + 2), T®_ hn commutes with on DFS„(A + 2). Therefore the commutator of T^_ lTl 

with elements in H^r^ annihilates all states not in DFS n (AT).[] To show that commuting T^-i n with H^T^ generates 
suitx) on the T-states of DFS„(A) we need the following lemma, which shows how to form su(d) from an overlapping 
su(d— 1) and su(2): 

Enlarging Lemma — Let TL be a Hilbert space of dimension d and let \i) £ TL. Assume we are given a set of 
Hamiltonians Hi that generates su(d— 1) on the subspace of TL that does not contain \i) and another set H2 that 
generates su(2) on the subspace of TL spanned by {[«), [,?)}, where \j) is another state in TL. Then [Hi,H2] (all 
commutators) generates su(d) on 7i under closure as a Lie-algebra (i.e., via scalar multiplication, addition and Lie- 
commutator). 

Proof — See Appendix |c] 

Now consider two states \i), \j) eDFS„(A) such that |i) ends in |11) and \j) is a T-state, but does not end in |11). 
Then we can generate su(2) on the subspace spanned by {\i), \j)} as follows: (i) We use the exchange interaction 
Xy = \i'){j'\ + [a prime indicates the bitstring with the last bit (a 1 in this case) dropped] in suitx) € H^T, x to 

generate a simultaneous action on DFS n (Af) and DFS„(AT + 2). This interaction is represented by a 2 x 2 cr^-matrix 
in the subspace spanned by {\i), \j)}. (ii) T^_ 1 n is represented by the 2x2 matrix diag(l, 0) = 5 (I + or z ) in the 
same subspace, and commutes with Xy on DFS„(A + 2) (since Xjj affects only B-states in DFS„(A + 2), and T^_ 1 n 
is non-zero only on states ending in |11)). Thus we can use it to create an independent action on DFS n (AT) alone: 
Yy = i[T n _ 1 ^,Xjj], Zjj = ^[YijjXjj]. 

Together {Yy, } generate su(2) independently on {\i), \j)} G DFS„(A). Since these operators vanish everywhere 
except on DFS„(A), their commutators with elements in [acting as su(tK)] will annihilate all other DFSs. 

Therefore, using the Enlarging Lemma, in this way all operations in su(tx) acting on DFS„(A) only can be generated. 

So far we have shown how to obtain an independent su(tx) on the T-states of DFS„(A) using H^r^ (for K < n— 4). 
To obtain an independent su{ok) on the B-states of DFS„(A) we use Hamiltonians in H^l x (acting on DFS n _i(A — 1) 
- the DFS from below). This will generate a simultaneous su(6a") in DFS„(A) and su{tK-2) in DFS„(A — 2). To 
eliminate the unwanted action on DFS n (AT— 2) we apply the previous arguments almost identically, except that now 



6 The argument thus far closely parallels the discussion above showing how to generate an independent su(2) on the 
{jOll), 1 101) } subspace of DFS 3 (-1), starting from the su(2) on DFS 2 (0) and T%. 
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we use T^_ 1 n to wipe out the action on all states except those ending in 1 00) . We thus get an independent su(ok) on 
DFS n (K). Together, the "above" and "below" constructions respectively provide independent su(i#-) and su(bfc) on 
DFS n (K). Finally, note that we did not really need both and T^, since once we established independent action 
on the T-states, we could have just subtracted out this action when considering the B-states. Also, the specific choice 
of T^ ,( ^ was rather arbitrary (though convenient): in fact almost any other diagonal interaction would do just as well. 

2. Mixing 

In order to induce operations between the two sets of paths (from "above" and from "below" ) that make up DFS„ (K) 
consider the effect of E„_i >n . This gate does not affect any paths that "ascend" two steps to (n,K) (corresponding 
to bitstrings ending in 1 00) ) and paths that "descend" two steps (ending in |11)), but it flips the paths that pass from 
(n - 2, K) via (n -l,K+l) with the paths from (n - 2, K) via (n -l,K-l) [see Fig. (|)]. It does this for all DFSs 
simultaneously. 

In order to get a full su(<Ik) on DFS n (K) we need to "mix" su(tK) (on the T-states) and su(6ff) (on the B-states) 
which we already have. We show how to obtain an independent su(2) between a T-state and a B-state. By the 
Enlarging Lemma this generates su(dx)- 

Since n > 3 DFS n (JC) contains states terminating in 1 00) and/or |11). Let us assume, w.l.o.g., that states termi- 
nating in |00) are present, and let \i) be such a state (B-state). Let \j) be a B-state not terminating in |00), and let 
|fc) = En— (\k) is a T-state). Let = — €= su{ok), and recall that we have independent su(bx)- 

Then as is easily checked, i[E„_i !n , Zy] = ~Yjk yields o~ y between \ j) and \k) only.Q In addition, Zj^ = ^[E n _i in , Yjfc] 
gives a z between |j) and |fc), thus completing a generating set for su(2) on the B-state |j) and the T-state |fc), that 
affects these two states only and annihilates all other states. This completes the proof. 

To summarize, we have shown constructively that it is possible to generate the entire Lie algebra su(dx) on a given 
weak collective-decoherence DFS„(if) of dimension dx, from the elementary composition of the operations of scalar 
multiplication, addition, Lie-commutators (conjugation by unitaries was not necessary in the WCD case). Moreover, 
this su(dx) can be generated independently on each DFS, implying that universal quantum computation can be 
performed inside each DFS„(i4T). Naturally, one would like to do this on the largest DFS. Since given the number 
of qubits n the dimensions of the DFSs are dx = ( K ) , the largest DFS is the decoherence-free subspace K — 0. In 
principle it is possible, by virtue of the independence result, to universally quantum compute in parallel on all DFSs. 

D. State Preparation and Measurement on the Weak Collective Decoherence DFS 

To make use of a DFS for encoding information in a quantum computer, in addition to the universal quantum 
computation described above, it must also be possible to initially prepare encoded states and to decode the quantum 
information on the DFS at the end of a computation. Encoding requires that the density matrix of the prepared 
states should have a large overlap with the DFS. Note that it is not necessary to prepare states that have support 
exclusively within the DFS. This follows from the fact that in our construction, while a computation is performed 
there is no mixing of states inside and outside of the DFS. If an initially prepared state is "contaminated" (has some 
support outside the DFS we want to compute on), then the result of the computation will have the same amount of 
contamination, i.e., the initial error does not spread. 

For example, suppose we can prepare the state p = (1 — p)\ip) {ip\ + p\tp±) (tp± | where is a state of a particular 
DFS and \ip±) is a state outside of this DFS. Then the computation will proceed independently on the DFS and the 
states outside of the DFS. Readout will then obtain the result of the computation with probability 1 — p. Repeated 
application of the quantum computation will give the desired result to arbitrary confidence level. 

There are many choices for the initial states of a computation and the decision as to which states to prepare should 
be guided by the available gates and measurements and the accuracy that is achievable. For efficient computation 
one should try to maximize the overlap of the prepared state with the desired initial DFS state. 

For the WCD case preparation of initial pure states is very simple. Suppose we are concerned with the S z error 
WCD-DFS. Pure state preparation into such a DFS then corresponds to the ability to prepare a state which has 
support over states with a specific number of |0) and |1) (eigenstates of the o~ z operator). This is particularly simple 
if measurements in the a z basis (|0) , |1)) as well as a x gates (to "flip" the bits) are available. 



7 Since E n _i jn = + \k}(j\ + \ + O, where O is some action on an orthogonal subspace. 
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The second crucial ingredient for computation on a DFS (in addition to preparation) is the decoding or readout 
of quantum information resulting from a computation. Once again, there are many options for how this can be 
performed. For example, in the WCD case one can make a measurement which distinguishes all of the DFSs and all of 
the states within this DFS by simply making a measurement in the a z basis on every qubit. Further, all measurements 
with a given number of distinct eigenvalues can be performed by first rotating the observable into one corresponding to 
a measurement in the computational basis (which, in turn, corresponds to a unitary operation on the DFS) and then 
performing the given measurement in the cr z basis, and finally rotating back. There are other situations where one 
would like to, say, make a measurement of an observable over the DFS which has only two different eigenvalues. This 
type of measurement can be most easily performed by a concatenated measurement [32] ]. In this scheme, one attaches 
another DFS to the original DFS, forming a single larger DFS. Then, assuming universal quantum computation over 
this larger DFS one can always perform operations which allow a measurement of the first DFS by entangling it with 
the second DFS, and reading out (destructively as described for the WCD above) the second DFS. For example, 
suppose the first DFS encodes two bits of quantum information, |fc, 1}l, k, I — {0, 1}, and the second DFS encodes a 
single bit of quantum information {|0)l, Then one can make a measurement of the observable a z (8 I on the 

first DFS by performing an encoded controlled-NOT operation between the first and the second DFS, and reading 
out the second DFS in the encoded a z basis. For the WCD case the ability to make this destructive measurement on 
the ancilla (not on the code) simply corresponds to the ability to measure single a z operations. 

Finally, we note that for a WCD-DFS there is a destructive measurement which distinguishes between different 
DFSs (corresponding to a measurement of the number of |l)'s). One can fault-tolerantly prepare a WCD-DFS state 
by repeatedly performing such a measurement to guarantee that the state is in the proper DFS. The concatenated 
measurement procedures described above for any DFS are naturally fault-tolerant in the sense that they can be 
repeated and are non-destructive |3^,|]|. Thus fault-tolerant preparation and decoding is available for the WCD- 
DFS. 

VII. STRONG COLLECTIVE DECOHERENCE 

Strong collective decoherence on n qubits is characterized by the three system operators S^, S y and S z . These 
operators form a representation of the scmisimple Lie algebra su(2). The algebra A generated by these operators can 
be decomposed as0 

n/2 

A= I nj ®gl{2J+l,€) (7.1) 

J=0(l/2) 

where J labels the total angular momentum of the corresponding Hilbert space decomposition (and hence the or 
1/2 depending on whether n is even or odd respectively) and gl(2J + 1, C) is the general linear algebra acting on a 
space of size 2 J + 1. The resulting decomposition of the system Hilbert space 

n/2 

H s = C nj <g> C 2l/+ i (7.2) 

J=0(l/2) 



is exactly the reduction of the Hilbert states into different Dicke states |32 53 . The degeneracy for each J is given by 
©: 

= (2J+l)n! 
J (n/2 + J + l)!(n/2 — J)! ' 1 ' ' 



Eq. (7.1) shows that given J, a state | J, A, /i) is acted upon as identity on its A component. Thus a DFS is defined by 
fixing J and /i. As we will show later, A corresponds to the paths leading to a given point (n, J) on the diagram of 
Fig. (|). 

The DFSs corresponding to the different J values for a given n can be computed using standard methods for the 
addition of angular momentum. We use the convention that |1) represents a \j = 1/2, mj = 1/2) particle and |0) 



Note that as a complex algebra {S x , S y ,S z } span all of gl(2), not just su(2). 
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represents a \j = 1/2, rrij = —1/2) particle in this decomposition although, of course, one should be careful to treat 
this labeling as strictly symbolic and not related to the physical angular momentum of the particles. 

The smallest n which supports a DFS and encodes at least a qubit of information is n = 3 |3^] . For n — 3 there are 
two possible values of the total angular momentum: J = 3/2 or J = 1/2. The four J = 3/2 states | J, A, fi) = |3/2, 0, //) 
(/i = in j = ±3/2, ±1/2) are singly degenerate; the J = 1/2 states have degeneracy 2. They can be constructed by 
either adding a J%2 = 1 (triplet) or a J12 = (singlet) state to a J3 = 1/2 state. These two possible methods of 
adding the angular momentum to obtain a, J — 1/2 state are exactly the degeneracy of the algebra. The four J = 1/2 
states are: 



1 1,0,0) 

lio,i) 



|0,0) ® 
|0,0) <g> 



II 1 

I 2 ' 2 



1 



(|010) -|100)) 



^(|011) -|101)) 



Hi) = 



11,1,0) = 

lii,i) = 



^(V2|l,l) 



)®\h-h)-\ 



- 10,0) ® 
1,0) <g> 



1 

2 ' 2 



-i» = 

5» 



1 

^(2|110) 



(-2|001) -1 
101 



|010) + |100» 
-|011)) 



(7.4) 



where in the first column we indicated the grouping forming a logical qubit; in the second we used the |J, A, fJ,) 
notation; in the third we used tensor products of the form \ Ji2,rnj 12 ) ® l^wjg); and in the fourth the states are 
expanded in terms the single-particle \j = 1/2, rrij = ±1/2) ba sis us ing Clebsch-Gordan coefficients. These states 
form a decoherence-free subsystem: the decomposition of Eqs. (7.1), ( [772]) ensures that the states {|^, 0, 0), ||, 0, 1)} 
are acted upon identically, and so are the states {||, 1,0), ||, 1, 1)}. Thus information of a qubit a\0z) + /3|1) should 
be encoded into these states as 



[(a*|0 L )+/3*|l L ))(a(0 L |+/3(l L |)][ 7 oo|0)(0|+7oi|0)(l|7oi±7io|l)(0|+7ii|l)(l|] 



(7.5) 



where 7^ form the components of a valid density matrix (unity trace and positive). Using Eq. (7.1) It follows that each 
of the S Q 's act on p in such a manner that only the A component is changed. Indeed, the S Q 's act like a corresponding 
a a in the /i-basis because this basis is two-dimensional, and a a are the two dimensional irreduc ible representations of 
su(2). These considerations are illustrated in detail for the exchange interaction in Sec. VII C . 



The smallest decohcrencc-frec subspace (as opposed to subsystem) supporting a full encoded qubit comes about 
for n — 4. Subspaces for the SCD mechanism correspond to the degeneracy of the zero total angular momentum 
eigenstates (there are also two decoherence-free subsystems with degeneracy 1 and 3). This subspace is spanned by 
the states: 



|0 L ) = 10,0,0) = |0,0) (8 |0,0) = -|(|01) - |10))(|01) - |10)) 
|1 L ) = |0, 1,0) = -L(|l, 1) O |1,-1) - |1,0) ® |1,0) + |1, -1) <g) |1, 1)) 
(2|0011) + 2|1100) - |0101) - 1 1010) - |0110) - |1001>). 



(7.6) 



The notation is the same as in Eq. (7.4), except that in the second column we used the notation | J12, mj 12 )® | J34, mj 34 ) 
which makes it easy t o se e how the angular momentum is added. 

As seen from Eqs. (7.4) and (|7.6|), there is a variety of useful bases which one can choose for the SCD-DFSs. We 
now show how the generic basis |J, A,/x) can be given both a graphical and an angular momentum interpretation. 
Consider the addition of angular momentum as more particles are included, similar to the construction we used in 
the WCD case. To construct the n qubit SCD-DFS for a specific J, denoted in this section as DFS n (J), one takes 
DFS„_i(J— 1/2) and DFS„_i( J+l/2) and uses the angular momentum addition rules to add another qubit (j = 1/2). 
Table (Q) presents the degeneracy of the J th irreducible representation for n qubits. The entries are obtained just as 
in Pascal's triangle, except that half of the triangle [the bottom according to the scheme of Table (ffl)] is missing. 



j — s 

J — 2 

J = 2 

J — 3 

J — 2 

J = 1 



TABLE I. Strong collective decoherence DFS dimensions. 
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J= i 1 2 5 

J=0 1 2 5 

n = 1 n = 2 n = 3 n = 4 n = 5 n = 6 
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Table (Q) demonstrates how the degeneracies of the (n — l)-qubit J± 1/2 irreducible representations (irreps), i.e., 
the dimensions of DFS„_i(J ± 1/2), add to determine the dimension of DFS„(J). This method of addition of the 
angular momentum leads to a natural interpretation of the | J, A, fj,) basis for the SCD-DFSs which we now present. 

Define the partial collective operators 

i=l 

This can be used to find a set of mutually commuting operators for the SCD-DFSs: the partial total angular momentum 
operators 

(S fc ) 2 = £ (S^) 2 . (7.8) 

a— x,y ,z 

As shown in Appendix 

[(S fc ) 2 , (S 1 ) 8 ] -0 Vk,l. (7.9) 

Thus the {(S^) 2 } can be used to label the SCD-DFSs by their eigenvalues Jk- 

In order to make the connection between the addition of angular momentum and the Dicke states one should, 
however, use 

i=l 

With this definition (s fc ) 2 = X) q ( s q) 2 i s i us ^ the operator whose eigenvalue for the J th irrep of the k qubit case is 
Jk{Jk + 1)- We label the basis determined by the eigenvalues of (s fc ) 2 by 

\Ji, J2, J3, ■ ■ ■ , J n -i,J;mj), (7-11) 

where 

(s fc ) 2 | Ji, J 2 , J 3 , . . ■ , Jn-x, J; mj) = J k (Jk + 1)1 Ji, -h, J3, ■ ■ ■ , J n -i, J] mj), (7.12) 

and where for consistency with the \J, A, //) notation we use J for J„. As in the WCD case, the degeneracy which leads 
to the SCD-DFS can be put into a one-to-one correspondence with a graphical representation of the addition of angular 
momentum. Here, however, each step does not simply correspond to adding a |0) or |1) state but instead corresponds 



J — 1/2 1 particle (note the 
) the horizontal axis counts 
qubit [note the similarity 



to combining the previous spin J particle with a spin 1/2 particle to create a J + 1/2 or 
absolute value so that the total spin is positive). In the graphical representation of Fig. (n 
qubits, and the vertical axis corresponds to the total angular momentum Ji up to the v 

to Table (jlj)]. Each SCD-DFS state then corresponds to a path constructed by successively moving up or down 1/2 
unit of angular momentum, starting from a single qubit with J\ = 1/2 . For example, the two DFS3(l/2) states are 
{ 1 1/2, 0, 1/2; ±1/2), 1 1/2, 1, 1/2; ±1/2)} (corresponding, respectively, to the paths "up, down, up" and "up, up, down" 
and mj 3= i/2 = ±1/2), and the two DFS 4 (0) states are {|l/2, 0, 1/2, 0; 0), 1 1/2, 1, 1/2, 0; 0)}. Clearly, the set of paths 
J„ = { Ji, J2, J3, . . . , Jn-ii Jn} with fixed J„ counts the degeneracy of DFS n (</„). Therefore we can identify the 
general degeneracy index A (of \J,X,/i}) with J„. Similarly, the dimensionality index /1 can now be identified with 
mj n . Finally, as claimed above J is just the final J n . 



A. The Stabilizer and Error Correction Properties 

Note from Eq. (|7.7|) that the system operators S a = SJJ. Therefore they can only affect the last component \J n ; mj n ) 
of the DFS states. By the identification of the degeneracy index A with the paths { Ji, . . . , J n -i, J}, and from the 
general expression ( |7.l|) for the action of the S Q , we know that S a acts only on the dimensionality component: 

S a | Ji, . . . , J n -i,J;mj) = I Ji, . . . , J„_i, J) (g) (P a \mj)), (7-13) 

where the P Q are a 2 J + 1 dimensional representation of su(2) acting directly on the \mj) components of the DFS. 
The corresponding DFS stabilizer is 
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T)(v) = T)(v x ,v y ,v z ) = exp 

_a— x,y : z 

For the J = DFSs this reduces to all collective rotations+contractions 



v a {S a -I®P a ) 



(7.14) 



D(fl) 



exp 



,a=x,y 1 z 



) exp [v ■ <Ji 



I cos I 



■ sin \v\ 



where ||u|| = (X^a^a) 1 ^ 2 ma y be complex. Thus DFS„(0) protects against all processes described by Kraus oper- 
ators that are linear combinations of collective rotations+contractions exp [v ■ a] . The situation for J + 1 is more 
complicated to calculate analytically. 

Let us now comment briefly on the error-correction and detection properties of DFS„(0): The stabilizer elements 
are tensor products of identical 1-qubit operators, including the following elements of the Pauli group: X® ra , Y®" and 
Z®" . Thus, for any odd-multiple 2k — 1 < n of single qubit errors X, Y and Z there is an element in the stabilizer 
that anticommutes with it: The code can detect any such error. The J = SCD-DFS is an error correcting code of 
distance 2. 



B. Nontrivial Operations 



Are there any single-qubit operators which preserve a SCD-DFS (and thus allow for nontrivial operations on the 
DFS)? There are no nontrivial single-qubit operators that commute with all S Q operators, since 

= i/^ijEgp-yol, (7-15) 
i i 

which vanishes iff a = /3. Therefore there are no single-qubit operators which preserve all SCD-DFSs simultaneously. 

As for two-qubit operators, the only such Hermitian operators which commute with the S a are those that are 
proportional to the exchange interaction [Eq. ( |6.9| )]: Eij\k)i\l)j = \l)%\k)j, where i,j label the qubits acted upon [ p9[ . 
In both the single- and two-qubit cases, there could be additional operators in the generalized commutant T (e.g., 
for n = 4 qubits there is an operator which mixes the different J's and preserves DFS^O): T = | J = X, Ai, Hx){J = 
2, Ai,/ii| + H.c). We will not be concerned with such operations as they are not needed in order to demonstrate 
universality, and since we will show that the exchange operator is sufficient for any SCD-DFS. Our task is thus to 
show that exchange interactions alone suffice to generate the entire SU(N) group on each iV-dimensional DFS, in the 
SCD case. 



C. Quantum Computation on the n — 3 and n — 4 qubit SCD-DFS 



We begin our discussion of universal quantum computation on SCD-DFSs by examining the simplest SCD-DFS 
whic h su pports encoding of quantum information: the n = 3 decoherence-free subsystem. We label these states as in 
Eq. (7.4) by \ J, A,/x). Recall that the J = 3/2 irrep is not degenerate and the J = 1/2 irrep has degeneracy 2. The 
J = 3/2 states can be written as ||, 0, fx), with /i = mj = ±3/2, ±1/2. Since the action of exchange does not depend 
on fi (recall that it affects paths, i.e., the A component only) it suffices to consider the action on the representative 
/i = 3/2 only: 1 1 1 1 ) . Let us then explicitly c alcu late the action of exchanging the first two physical qubits on this 
state and the four J = 1/2 states. Using Eq. (7.4): 



Ei 2 |-,0, -) 
1 2 2 7 

1 



Ei 2 |-,0,1) 



E ia |lll) = ||,0, |) 
1 



E 12 |-,0,0) = E 12 — (|010) - |100» = — (|100) - |010)) = -|-,0,0) 



V2 



V2 



E 12 -L (|011) - |101)) = ± (|101) - |011)) = 



1 



1 



E 12 |-, 1, 0) = E 12 ^= (-2|001) + |010) + |100)) = |-, 1, 0) 
E 12 |~ 1, 1) - E 12 -L (21110) - |101) - |011)) = |i, 1, 1). 



(7.16) 
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Focusing just on the J — 1/2 states, the exchange action on |A) (g) \n) can thus be written as: 



E12 = -o z ® I. (7.17) 

Since the action of the S a operators on the J = 1/2 states is I ni , 2 <8> gl(2) according to Eq. (7.1), this explicit form 
for E12 confirms that is has the expected structure of operators in the commutant of the algebra spanned by the S a . 
It can also be seen that quantum information should be encoded in the |A) component, as discussed before Eq. ([7.5| ). 

Using similar algebra it is straightforward to verify that the effect of the three possible exchanges on the n = 3 
DFS states is given by: 














El2 — 





-1 
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E23 — 
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1 
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j E13 = 





1 

2 


V3 




Vo 


2 


2 / 



(7.18) 




where the rows and columns of these matrices are labelled by the basis elements {|J = 3/2, A = 0),|J = 1/2, A = 
0), I J = 1/2, A = 1)}. As expected from general properties of the commutant, the exchange operators do not mix the 
different J irreps. Now, 

-(E12 + E13 + E23) = 
-(— E12 + E13 + E23) = 

-L(E 13 -E 23 ) = I 1 I , (7.19) 

showing that the last two linear combinations of exchanges look like the Pauli a z and a x on DFS 3 (l/2). Using a 
standard Euler angle construction it is thus possible to perform any SU(2) gate on this DFS. Moreover, it is possible 
to act independently on DFS3(3/2) and DFS3(l/2). In other words, we can perform U(l) on DFS3(3/2) alone, and 
SU(2) on DFS 3 (l/2) alone. Note, however, that at this point we cannot yet claim universal quantum computation 
on a register composed of clusters of DFS3( J)'s (J constant) because we have not shown how to couple such clusters. 

For n ~ 4 the Hilbert space splits up into one J = 2-irrep [DFS4(2)], three J = 1-irreps [DFS^l)], and two 
J = 0-irreps [DFS4(0)] - see Tabic (||). Direct calculation of the effect of exchange on these DFSs shows that we can 
independently perform su(l) (i.e. zero), su(3), and su(2). In particular, we find that |2^ , |32]| : 

1 % 2 

X = -^=(E 2 3 — Ei 3 ) Y = ^-^=[E 2 3 — E13, E 34 ] Z = -[Y, X] = — E12 (7.20) 

act as the corresponding su(2) Pauli operators on DFS4(0) only. Further, the following operators act independently 
on the J — 1-irreps (rows and columns are labelled by A — 0, 1,2. The action occurs simultaneously on all three fi 
components corresponding to a given A): 

„■ /0 - , , 

Yi3 = — =[Ei2,E34]= , X 13 = -[E 12) Y 
2 v 2 \i I 

i ( 1 

2 13 lw 





J[Y 13) X 13 ]= 10 I , Y 2 3 = ^[E23,Z 13 ]= I -i I . (7.21) 



2 

These operators clearly generate su(3), and hence we have an independent SU(3) action on DFS4(1). 

D. Universal Quantum Computation on the n > 5 qubit SCD-DFSs 

We are now ready to prove our central result: that using only the two-body exchange Hamiltonians every unitary 
operation can be performed on a SCD-DFS. More specifically: 
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Theorem 5 — For any n > 2 qubits undergoing strong collective decoherence, there exist sets of Hamiltonians H ™ 
obtained from exchange interactions only via scalar multiplication, addition, Lie-commutator and unitary conjugation, 
acting as su(dj) on the DFS corresponding to the eigenvalue J. Furthermore each set acts independently on this DFS 
only (i.e., with zeroes in the matrix representation corresponding to their action on the other DFSs). 

In preparation for the proof of this result let us note several useful facts: 

(i) The excha nge operators do not change the value of mj, because they are in the commutant of A. = {S a } [and 
recall Eq. ( 7.13| )]. Therefore in order to evaluate the action of the exchange operators on the different DFS n (J) ( 



n 



given) it is convenient to fix mj, and in particular to work in the basis given by the maximal mj value (mj = J). 
Expressions for these "maximal" states in terms of | J\, J%, . . . , J n -2 \ mj) and the single qubit states of the last two 
qubits are given in Appendix [b|. 

(ii) Every (s fe ) 2 can be written as a sum of exchange operators and the identity operation EjJ. This follows from 



Eq. (Al) and noting that the exchange operator can be expanded as 

1 



E« = g I 1 + « + °>l + «) > ( 7 - 22 ) 



so that: 



^f = k(l-^)l+ \ £ E , (7.23) 

Thus (s fe ) 2 is a Hamiltonian which is at our disposal. 

We are now ready to present our proof by induction. Recall the DFS-dimensionality formula for nj, Eq. fl7.3p. We 
assume that it is possible to perform su(nj) on each of the different DFS„_i(J) independently using only exchange 
operators and the identity Hamiltonian. Our construction above proves that this is true for 3 and 4 qubits. The 
assumption that the actions we can perform can be performed independently translates into the ability to construct 
Hamiltonians which annihilate all of the DFSs except a desired one on which they act as su(nj). 

As in the WCD case a specific DFS n (J) of dimension nj splits into states which are constructed by the subtraction of 
angular momentum from DFS n _i(J+l/2) (T-states), or by the addition of angular momentum to DFS„_i( J— 1/2) (In- 
states) [see Fig. (g)]. Performing su(n J+1 / 2 ) on DFS„_i(J+l/2) will simultaneously act on DFS„(J) and DFS„( J+l). 
In other words, sit(n,/ +1 / 2 ) on DFS„_i(J + 1/2) acts on both the B-states of DFS„(J + 1) and on the T-states of 
DFS„( J). We split the proof into three steps. In the first step we obtain an su(2) set of operators which acts only on 
DFS„(J) and mixes particular B- and T-states. In the second step we expand the set of operators which mix B- and 
T-states to cover all possible su(2) algebras between any two B- and T-states. Finally, in the third step we apply a 
Mixing Lemma which shows that we can obtain the full su(nj) (i.e., also mix B-states and mix T-states). 

1. T- and B-Mixing 

There are two simple instances where there is no need to show independent action in our proof: (i) The (upper) 
J = n/2 -irrep is always 1-dimcnsional, so the action on it is always trivial (i.e., the Hamiltonian vanishes and hence 
the action is independent by definition); (ii) For odd n the "lowest" DFS„(l/2) is acted upon independently by 
the su(jiq) from DFS n _i(0) [i.e., su(jiq) cannot act "downward"]. In order to facilitate our construction we extend 
the notion of T and B-states one step further in the construction of the DFS. TB-states are those states which are 
constructed from T-states on (n — l)-qubits and from the B-states on n-qubit states [see Fig. (|)]. Similarly we can 
define the BT, TT, and BB-states: 



TT) = 




■ 7 J 71 


|BT) = 




• ) Jfl 


TB) = 


\Ju- 


■ ; J 71 


|BB) = 


\Ju- 


■ 1 J 71 



(7.24) 



(7.25) 



Every DFS n (J) can be broken down into a direct sum of TT, BT, TB, and BB-states; e.g., as seen in Fig. (Q), in 
DFS 6 (1) there are 1 TT, 3 TB, 3 BT and 2 BB states. Note that for J = n/2 — 1 there are no TT-states, for J = 



2G 



there are no BB and BT-states, for J —1/2 there are no BB-states, and otherwise there are as many TB as there are 
BT states, 



At this point it is useful to explicitly give the action of exchange on the last two qubits of a SCD-DFS. Using Eq. (BS ) 
we find (assuming the existence of the given states, i.e., n large enough and J not too large) the representation 



E 



n ,n—l — 



1 \ TT 

_ | -cos(# J+1 ) sm(9 J+1 ) BT 

sin(0 J+ i) cos((9 J+ i) TB 

1 / BB 



(7.26) 



where tan(#j) = 2y/ J(J + 1). Thus exchange acts to transform the BT and TB states entering a given DFS into 
linear combinations of one another, while leaving invariant the BB and TT states. 

Let us now consider the action of su(nj_ 1 / 2 ) from DFS„_i(J — 1/2) [see Fig. (|J)]. It acts on DFS n (J — 1) and 
DFS„(J) simultaneously. However, since the T-states of DFS„(J — 1) and the B-states of DFS„(J) share the same 
set of quantum numbers {Ji, J n _i}, the action of the su(nj— 1/2) operators is identical on these two sets of states. 

We first deal with the case where the number of BT-states of DFS„(J) is greater than 1. As can be inferred from 
Fig. (Q), this condition corresponds to J < rt/2 — 1 and n > 4. We will separately deal with the J = n/2 — 1 case at 
the end of the proof. Let \a) and \b) be any two orthogonal BT-states of DFS„( J) (i.e., states differing only by the 
paths on the first n — 2 qubits). Corresponding to these are {\a'), \b')}: a pair of orthogonal BT-states of DFS„(J). 
One of the elements in su(nj_i/ 2 ) is the traceless operator C = \a)(a\ — \b){b\, which we have at our disposal by the 
induction hypothesis. Consider i[E ni „_i, C]: since E„.„_i acts as identity on BB states, even though C has an action 
on DFS„( J — 1) the commutator acting on the BB states of DFS„ ( J — 1) vanishes. The action of i[E njn _i, C] on the 
BT and TB states can be calculated by observing, using Eq. ( 7.26 ), that the matrix representations of C and E nitt _i 



are, in the ordered {\a'}, \b'\ 



,\b}} basis: 



C = diag(0, 0,1,-1) = -(I ( 



E 



n.n—1 



cr z - a z 
sin(0j) 



f - cos(Oj) sm(e.j) 

-cos(0j) sin(0j) 

sin(6»j) cos(6»j) 

V sin(e». 7 ) cos(0j) 



= — cos(0j)cr z ® I + sm(8j)cr x ® I. 



This yields: 



i S m(ej)(-\a)(a'\ + \a')(a\ + \b)(b'\-\b'}(b\). 



(7.27) 



(7.28) 



Now let \c) be a TT-state of DFS„( J). Such a state always exists unless J = n/2 — 1, which is covered at the end of 
the proof. Then there is an operator D = |a')(a'| — \c)(c\ in su(n j+1/2)® It follows that: 



X„ 



1 



-i[i[B n , n _i, C],D] = \a)(a'\ + \a')(a\ 



(7.29) 



sin(0j) 

acts like an encoded a x on \a) and \a!) and annihilates all other states. Further, one can implement the commutator 

Y 0a / = i[X oo ,,D]=i(|a)<a / |-|o / )(o|) ) (7.30) 

which acts like an encoded a y on |a) and \a'). Finally, one can construct Z QQ ' = i[X aa /, Y 00 /] = |a)(a| — |a')(a'|. Thus 
we have shown that for J < n/2 — 1 we can validly (using only exchange Hamiltonians) perform su{2) operations 
between \a), a specific B-state and \a'), its corresponding T-state, on DFS„(J) only. 



2. Extending the su(2) 's 



We now show that by using the ope rati on of conjugation by a unitary we can construct su(2) between any two 
B and T-states. To see this recall Eq. (4.6), which allows one to take a Hamiltonian H and turn it via conjugation 



We need to subtract |c){c| in order to obtain a traceless operator. 
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by a unitary gate into the new Hamiltonian H c ff = UHU*. By the induction hypothesis we have at our disposal 
every SU gate which acts on the T-states of DFS„(J) [and simultaneously acts on the B-states of DFS„(J + 1)] and 
also every SU gate which acts on the B-states of DFS„( J) [and simultaneously acts on the T-states of DFS„( J — 1)]. 
Above we have shown how to construct X, Y, and Z operators between specific T- and B-states: |o') and \a). 
Let \i') and \i) be some other T- and B-states of DFS„(J), respectively. Then we have at our disposal the gate 
P»'i = I a') + 10 + |o)(i| + |i)(o| + O where O is an operator which acts on a DFS other than DFS ra ( J) (included 
to make Pj/, an SU operator). It is simple to verify that 

= P^X^Pj,; = + (7.31) 

which acts as an encoded a x between \i') and \i). Note that because X Qa / only acts on DFS„(J), Xj^ will also only 
act on the same DFS. Similarly one can construct Yj', = Pi'tY ao /Pj^ and Zj/j = P^Z^/P^ which act, respectively, 
as encoded o~ y and o~ z on and Thus we have shown that one can implement every su(2) between any two T- 
and B-states in DFS n (J). Each of these su(2) operations is performed independently on DFS„(J). 



3. Mixing T- and B-States 

Next we use a Lemma proved in Appendix 

Mixing Lemma: Given is a Hilbert space TL — TL\ ® TL2 where AvcaTLj = nj. Let and {|«2)} be orthonormal 

bases for Tii and H.2 respectively. If one can implement the operators Xi 1 i 2 = |«i)(?2| + |*2)(*i|, = i|ii)(i2| — 

«|i2)(ii|, and Z^i 2 = — 1*2) 1 , then one can implement su(n\ + n^) on H. 

Above we have explicitly shown that we can obtain every X^^, Y^^, and Z^^ acting independently on DFS„(J). 
Thus direct application of the Mixing Lemma tells us that we can perform su(nj) independently on this DFS. 

Special case of J = n/2 — 1: We have neglected DFS„(n/2 — 1) because it did not contain two different BT-states 
(nor a TT) state. The dimension of this DFS is n—1. We now show how to perform su(n — 1) on this DFS using the 
fact that we have already established su(nj =n /2-2) on DFS„(n/2 — 2). First, note that by the induction hypothesis 
we can perform su(nj =n / 2 -3/2) independently on DFS n _i(n/2 — 3/2). As above, this action simultaneously affects 
DFS„(n/2 — 1) and DFS„(n/2 — 2). However, since we can perform sit(nj = „/ 2 -2) on DFS„(n/2 — 2), we can subtract 
out the action of su(nj =n / 2 -3/2) on DFS„(n/2 — 2). Thus we can obtain su(nj =n / 2 -3/2) on all of the B-states of 
DFS n (n/2 — 1). But the exchange operator E ni „_i acts to mix the B-states with the single T-state of DFS n (n/2 — 1). 
Thus we can construct an su(2) algebra between that single-T state and a single B-state in a manner directly analogous 
to the above proof for J < n/2 — 1. Finally, by the Enlarging Lemma it follows that we can obtain su{n — 1) on 
DFS n (ra/2 - 1). 

This concludes the proof that the exchange interaction is universal independently on each of the different strong- 
collcctive-decoherence DFSs. 



E. State Preparation and Measurement on the Strong Collective Decoherence DFS 

At first glance it might seem difficult to prepare pure states of a SCD-DFS, because these states are nontrivially 
entangled. However, it is easy to see that every DF subspace contains a state which is a tensor product of singlet 
states: 

|0 o )=(^=J 8^(101) -|10)), (7.32) 

because these states have zero total angular momentum. Thus a supply of singlet states is sufficient to prepare DF 
subspace states. Further, DF subsystems always contain a state which is a tensor product of a DF subspace and a 
pure state of the form |1) ® ■ ■ ■ ® |1). This can be seen from Fig. (|j), where the lowest path leading to a specific 
DFS„(J) is composed of a segment passing through a DF subspace (and is thus of the form |0d)), and a segment 
going straight up from there to DFS„(J). The corresponding state is equivalent to adding a spin-0 (DF subspace) 
and a spin-J DF subsystem (the | J, mj = J) state of the latter is seen to be made up entirely of |1) ® • • • ® |1)). In 
general, addition of a spin-0 DFS and a spin-J DFS simply corresponds to tensoring the two states. Note, however, 
that addition of two arbitrary DF subsystems into a larger DFS is not nearly as simple: concatenation of two J / 
DFSs does not correspond to tensoring. 

Pure state preparation for a SCD-DFS can thus be as simple as the ability to produce singlet states and |1) states (it 
is also possible to use the | J, mj = — J) = |0) ® ■ ■ ■ <B> |0) or any of the other \ J, mj) states plus singlets). Other, more 
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complicated pure state preparation procedures are also conceivable, and the decision as to which procedure to use is 
clearly determined by the available resources to manipulate quantum states. The pure state preparation of singlets 
and computational basis states has the distinct advantage that verification of these states should be experimentally 
achievable. Such verification is necessary for fault-tolerant preparation [j43| . 

Measurements on the SCD -DFS can be performed by using the concatenated measurement scheme detailed in the 
WCD-DFS discussion [Sec. ( VI D| )l . In particular, by attaching a SCD-DF subspace ancilla via concatenation, one 
can construct any concatenated measurement scenario. All that remains to be shown is how to perform a destructive 
measurement on such an ancilla. In J|^,Q such schemes are presented for the n = 4 SCD-DF subspace which 
encodes a single qubit of information. We will not repeat the details of these schemes here, but note that they involve 
measurements of single physical-qubit observables and thus are experimentally very reasonable. Further, we note 
that the ability to perform a concatenated measurement scenario by concatenating an ancilla DFS composed of a 
single encoded-qubit, can be used to perform any possible concatenated DFS measurement scenario. As mentioned 
in the WCD case, the concatenated measurement procedures are fault-tolerant. Thus we have shown how to perform 
fault-tolerant preparation and decoding on the SCD-DFS. 



VIII. UNIVERSAL FAULT-TOLERANT COMPUTATION ON CONCATENATED CODES 

So far we have shown how to implement universal computation with local Hamiltonians on a DFS corresponding 
to a single block of qubits. This construction assumes that the only errors are collective. This is a very stringent 
symmetry requirement, which obviously becomes less realistic the larger the number of particles n is. It is thus 
desirable to be able to deal with perturbations that break the collect ive-decoherence (permutation) symmetry. To 
this end we have previously studied the effect of symmetry- breaking perturbations on decoherence-free subspaces [^ij , 
and have proposed a concatenation method to make DFSs robust in the presence of such perturbations. This method 
embeds DFS blocks of four particles (each block constituting a single encoded qubit) into a QECC p^] . The QECC 
in the outer layer then takes care of any single encoded-qubit errors on each of its constituent DFS-blocks; in fact the 
code can correct for any "leakage" error taking a state outside of the DFS, by transforming this into a single encoded 
qubit error on the QECC. By choosing an appropriate QECC it is thus possible to deal with any type of non-collective 
error on the encoded DFS-qubits. In particular, by using the "perfect" 5-qubit code it is possible to correct all 
independent errors between blocks of four particles. 

The problem with this construction so far was, that in order to correct on the outer QECC, it is necessary to 
perform encoded operations on the constituent DFSs in a fault-tolerant way, using (realistic) local interactions. 
Specifically, it is necessary to be able to implement all single encoded-qubit operations on the DFS-qubits of the 
outer QECC, as well as operations between two DFS-blocks (see for details). Given that one can perform single 
qubit (or "qupit" for higher-dimensional DFSs) operations on each DFS-block, the only additional gate necessary to 
implement error-correction and universal quantum computation on a concatenated QECC-DFS, is any non-separable 
two-encoded-qubit gate K between any four states in the two DFS-blocks. For instance, a controlled-phase operation, 
which gives a phase of —1 to \0l0l) and leaves all other states unchanged. In fact, it is sufficient to be able to perform 
this gate K between neighboring blocks only []5(| . 

The above results give us the tools to perform single DFS-qubit (or qupit) operations on a block. To construct 
an encoded K between two neighboring blocks, we assume that the corresponding physical qubits are spatially close 
together during the switching time of the gate. Since the symmetry of collective decoherence arises from the spatial 
correlation of the decoherence process, we can further assume that during this switching time, both DFS-blocks couple 
to the same bath mode. This assumption is physically motivated by the expectation that collective decoherence occurs 
in the analog of the Dicke limit of quantum optics, where the qubits have small spatial separations relative to the 
bath correlation length . Then the two DFS-clusters temporarily form a bigger DFS and we can use the universal 
operations we have constructed previously on this big DFS to implement the desired gate K. 

This makes the concatenated QECC-DFS fully workable as a code supporting universal fault-tolerant quantum 
computation. 



IX. SUMMARY AND CONCLUSIONS 



In this paper we have settled the issue of quantum computation with realistic (few-body) means on both decoherence- 
free subspaces and decoherence-free (noiseless) subsystems (DFSs) for two important forms of decoherence: collective 
phase damping ("weak collective decoherence"), and collective phase damping plus collective dissipation ("strong col- 
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lective decoherence" ) . This resolves an outstanding question as to whether universal computation on these physically 
relevant DFSs by using just 1- and 2-body Hamiltonians is possible. 

The implications of this result for the usefulness of DFSs are drastic. They put the theory of DFSs on an equal 
footing with the theory of quantum error correction, in that the full repertoire of universal fault tolerant quantum 
computation is now available on DFSs for collective decoherence: the most important pertinent decoherence process. 
Moreover, the strict assumption of collective decoherence can be lifted by allowing for perturbing independent qubit 
errors. As we proposed earlier it is possible to stabilize DFSs against such errors by concatenation with a quantum error 
correcting code (QECC). However, to be able to implement error-correction and fault-tolerant universal computation 
on these concatenated codes a crucial (and so far missing) ingredient was the ability to perform encoded operations 
on the DFS-blocks fault-tolerantly. This paper settles that matter, showing constructively that DFSs can be made 
robust. 

Furthermore, this paper reports on a general framework incorporating both DFSs and QECCs, and generalizes 
the theory of stabilizer codes to the (non-abclian) DFS-casc. This framework enabled us to identify the allowed 
operations on a DFS and to show that these operations can be performed while maintaining a very strong form of 
fault-tolerance: the states remain within the DFS during the entire switching time of the gate. Our formalism should 
be readily applicable for other non-additive codes. 

There is an interesting duality between QECCs which are designed to correct single (or greater) qubit errors and 
DFSs. In QECC the errors are all single body interactions. The QECC condition therefore implies that any one or 
two-body Hamiltonian must take codewords outside of the code space because these interactions themselves look like 
errors. QECCs must leave their codcspacc in order to perform quantum computation on the encoded operations. This 
means that QECCs must have gates which act much faster than the decoherence mechanism so that a perturbative 
treatment can be carried out. QECC can correct small errors but the price paid for this is that gates must be 
executed quickly (not to mention that fault-tolerant gates must also be used). DFSs on the other hand, do not have 
the requirement of correcting single qubit errors and we have found that a single two-body interaction (exchange) is 
sufficient to generate universal quantum computation fault-tolerantly. DFSs have larger errors but this allows for an 
economy of Hamiltonians. 

As corollaries to our results on weak and strong collective decoherence two additional properties of the corresponding 
DFS encodings appear: 

• One can work on all DFSs in parallel: Since we are able to implement SU(d n ) on each DFS n (n=number 
of particles) independently, we can in principle work on all DFSs in parallel. This means that we can encode 
quantum information into each of the DFSs and perform calculations (possibly different) on all of them at once. 

• For the strong collective decoherence case the exchange gate is asymptotically universal: It is well 
known that the encoding efficiency of the singlet space of the strong collective decoherence-DFS for large n 
approaches unity More precisely, let k be the number of encoded qubits in the singlet (J = 0) sector of a 
Hilbert space of n qubits, then 

k 3 log, n 

lim - = 1 - -— 22— . (9.1) 

rwoo n 2 n 



We have established that the exchange gate alone (with an irrational phase) implements universal computation 
on each DFS and on the singlet space in particular. Thus, we find that, for large n, in order to achieve universal 
computation with nearly perfect efficiency, all we need to be able to perform is the exchange interaction. This 
result is very promising from an experimental point of view, since the exchange interaction is prevalent whenever 
there is a Heisenberg coupling between systems j2^]3^]. We emphasize that regardless of the decoherence 
mechanism, this implies that universal quantum computation can be achieved "asymptotically" using a single 
gate HiJ]. We conjecture that there are many more such two-body interactions which similarly provide such 
"asymptotic universality" on their own. 
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APPENDIX A: THE PARTIAL COLLECTIVE ANGULAR MOMENTUM OPERATORS ARE A SET OF 

COMMUTING OBSERVABLES 



We prove here that the partial collective operators = Sq ' 2 ' "''^ = J2i=i a a form a commuting set and hence, a 
good operator basis. Note first that 



( sfc ) 2 = E E «■ 

a—x,y,z 



(Al) 



Thus 



[(S fc ) 2 ,(S') 2 i 



E E <<, E E « 

i.j — l a—x.y,z m,n— 1 {3=x,y,z 



(A2) 



Terms with a = (3 obviously commute. Further, terms with (m = i,n = j), (to = j, n = i), or (i ^ m, n, j ^ m,n), 
commute, so we need only consider (i = m,j ^ n), (i = n,j ^ to) or (i ^ m,j = n), (i ^ n, j = in). In addition, 
assuming w.l.o.g. that I > k, terms with m 1 n> k also commute. Thus we are left with 



[(s fe ) 2 ,(s ; ) 2 ] = 2 £ £ k>^w]+ 2 E E E [«>^; 

i,j=l n(^j) = l (i^a=x,y,z ij=l m(=£i) = \ f3^a=x,y,z 

Using the fact that [a l a a J a , a l & af\ = i^e^^^ and [a l a a J a , er£ ! a 3 ^ = i ^ 7 e^a^u^u^: 

[(S fe ) 2 ,(S') 2 ]=2^ 2 ^ W ;« + 2 E E E wi<#* 

*iJ=l n(^j) = l a,f3,-y={x,y,z} i,j = l m(^i) = l a,fl,y={x,y,z} 

and both sums vanish due to the antisymmetric property of £ a/ 3 7 . 



(A3) 



APPENDIX B: MAXIMAL-mj STATES OF THE STRONG COLLECTIVE DECOHERENCE DFS 

We show how to recursively express the n-particle total spin-J states in terms of (n — l)-particle states. Let us 
focus on DFS„(J) and in particular on the maximal-mj state in it: 



\ip) = \ Ji,...,J n -i,J;mj = J). 



(Bl) 



In general (J ^ 0,n/2) there are two kinds of states: bottom (\iP}b) and top (|V')t) ones. The angular momentum 
addition rule that must be satisfied for adding a single spin-i particle is that 



The B-state comes from adding a particle to the maximal toj state in DFS„_i(J — 1/2), which is: 

|B) = \ Ji,..., J„_ 2 , J-^mj,,., = J - i). 



(B2) 



There is only one way to go from |B) to \i/)}b, namely to add 1/2 to mj n _ 1 = J — \ in order to obtain to,/ = J. Thus 



>B = |B)|i,i), 



(B3) 



where |i, i) is the single-particle spin-up state. The situation is different for the T-state, which is constructed by 



I 2 ' 2 

adding a particle to 



|T±) = |Ji,...,J n _ 2 ,J+^;mj B _ 1 = ^±^>- 



(B4) 
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These two possibilities give: 

|^> T = «1T + >||, -|> + /3|T_>||, |>. (B5) 

To find the coefficients a and (3 we use the collective raising operator s + = s x + is y , where we recall that Sq = 
i Ylj—i cr l a . Since \ip) is a maximal-mj state it is annihilated by s + = s^* . Similarly, |T+) is annihilated by s^™ . 
Therefore, since s + = s+ ~^ + 

S+ |T + >||,-|> = |T + )|i,i) 
s+ |T_>|ii> = V2J+T|T + )|i I), 

where in the second line we used the elementary raising oper ato r formula 3+\j, m) = + 1) — m(m + 1)] \j, m+1) 
with j = J + | and m = J — |. Application of s + to Eq. (B5) thus yields: 



a + V2J+T/3 = (B6) 
Hence, up to an arbitrary phase choice, we find that 



The special cases of J = 0, n/2 differ only in that the corresponding DFSs support just T- and B-states, respectively. 
The calculation of the coefficients, therefore, remains the same. 

In a similar manner one can carry the calculation one particle deeper. Doing this we find for the maximal-mj states 
(provided they exist): 



|TT) = \J x ,...,J n ^J+l,J+^,J;mj = J) = J 1; . . . , J„_ 3 , J + 1; TOj „_ 2 = J + -I) 



' 2J+1 , (1 11 _1 1 _1 1 1, 

(2J + 2)(2J + 3)' n " 3 ' + ;mj - 2 ~ } ^2' 2 }l 2' 2 ) + l 2' 2 }l 2'2' 



2 1111 

(2J + 2)(2J + 3) |Jl ' ' ■ ' ' J "~ 3 ' J + 1;mj "" 2 = J 1)1 2' 2 }l 2' 2 } 



|BT) = |Ji,..., J„_ 3 , J, J+ J;mj = J) = -J ^y^l J i>- ■ • ,Jn-3,J;m Jn _ 2 = J)|i, ~}|^,-^) 

|Ji,..., J„_ 3 , J;m Jn _ 2 = J)^,-- 1 ' 



X /(2J + 2)(2J+1)' ' ' ' "2' 2"2'2' 

IJi,.--, J„-3, J;mj„_ 2 = J-l)\l;,7;)\~, h 



(2J+l)(2J + 2)' A ' ' "~ J ' ' J - 2 ,l 2'2 ,l 2'2 , 



|TB) = |Ji,..., J„_ 3 , J, J- ^ J;mj = J) = -y 2 J + 1 1 Jli • - ' ' J "- 3 ' J ; mj "~ 2 = J ^ ,_ ^^'^ 

1 . v .l 1. 

Ji,..., J„_3, J;m Jn _ a = J — 1) -, -) -, -) 



|BB) = |J 1 ,...,J w _ 3 ,J-l,J-i J;mj = J) = |Ji,...,J n _ 3 , J-l;m Jn _ 2 = J - ±). (B8) 

Caution must be exercised in using these expressions near the boundary of Tabic (||) , where some of the states may 
not exist. 
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APPENDIX C: PROOFS OF THE LEMMAS 



Enlarging Lemma — Let TL be a Hilbert space of dimension d and let \i) G TL. Assume we are given a set of 
Hamiltonians Hi that generates su(d — 1) on the subspace of TL that does not contain \i), and another set H2 that 
generates su(2) on the subspace of TL spanned by {\i), \j}}, where \j) is another state in TL. Then [Hi,H2] (all 
commutators) generates su{d) on TL under closure as a Lie-algebra. 

Proof. We explicitly construct the Lie-algebra su{d) with the given Hamiltonians. Let TL C TL be the d — 1 
dimensional subspace Hi acts on. Let us show that we can generate su(2) between \k) S TL and \i). 

Let X y ee + e H 2 and X ]k ee \j)(k\ + \k)(j\ e Hi. Then 

YifcSttXj-fc.Xij^-iliXfcl+ilfcX*! (CI) 

acts as cr j, on the states \i), \ k). Similarly 

X ik =i[Y ij ,X jk ] = \i)(k\ + \k)(i\ (C2) 

yields a x on the space spanned by \i), \k). These two operations generate su(2) on \i), \k) for all \k) in the subspace 
of TL that does not contain Now the Mixing Lemma gives the desired result together with the observation that 
there we only use elements in [Hi, H2]. 

Mixing Lemma — Consider the division of an n dimensional Hilbert space TL into a direct sum of two subspaces 
TL\ © TL2 of dimensions n\ and n 2 respectively. Suppose that \i n ) is an orthonormal basis for TL n . Then the Lie 
algebras generated by X ilii2 = |ii)(i 2 | + |«2)(«i|, ^i u i 2 = i\h){ii\ - i\fa){h\i and Z ilji2 = \ii){h\ - \h){i2\ generate 
su(n). 

Proof. We explicitly construct the elements of su(n). Consider i[Xj li j 2 , 'Yj 1 J2 ]. Clearly, if i\ 7^ 12 =/= ji =/= 32 this 
equals zero and if i\ = j% and i 2 = j% then this commutator is — Z-j 1( j 2 . If, however, i\ = j\ and i 2 7^ j 2 this becomes 

*[X il)<a ,Y ilA ]=-|< a )(j a |-|i 2 >(< 2 |. (C3) 

Similarly: 

»[X <1 , ial Y il , ia ] = |< 1 )0i| + |j 1 )(*i|. (C4) 
Thus every \ik){ji\ + \ji){ik\ is in the Lie algebra. Similarly, z[Xj li j 2 , 'X.j 1 J2 ] yields 

i[X n , l2 ,X nj2 ] = i|i 2 )0'2| - i\h){h\ 

i[Xi u i 2 ,X jui2 ] = i\h)(ji\ - i\h)(h\ (C5) 

Thus every i\ik){ji \ — i\ji){ik \ is in the Lie algebra. Taking the commutator of these with the \ik)(ji \ + \ji) fal operators 
finally yields every |ifc)(ji| — \ji)(ik\- Since su(n) can be decomposed into a sum of overlapping su(2)'s p% , the Lie 
algebra is the entire su(n), as claimed. 
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dim = 




B 



FIG. 1. Weak collective decoherence. The horizontal axis is the number of qubits; the vertical axis is the number of O's minus 
the number of l's in each state (J„). Each state in the standard basis thus corresponds to a path from the origin following the 
indicated arrows. The dimension of a DFS corresponds to the multiple pathways through which one can arrive at the same J„. 
The insert shows the matrix structure of operators acting on DFSs(3), in terms of Top (T) and Bottom (B) states. Note that 
there is only one T-state entering DFSs(3), whence the action of exchange is represented by a 1 x 1 block. 
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FIG. 2. Action of exchange is to simultaneously flip different paths. Matrices are the representations of E34 on DFS4 
(lower) and DFS 4 (2) (upper). 
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su{d K+x ): 





FIG. 3. The structure of the pathways connecting adjacent DFSs in the weak collective decoherence case. The action of the 
different su Lie algebras is indicated by heavy arrows. 
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dim = 




FIG. 4. Strong collective decoherence. The horizontal axis is the number of qubits and the vertical axis is the total angular 
momentum. Each state in the DFS is represented by a pathway from the origin along the arrows as indicated. The insert shows 
the matrix structure of operators acting on DFSe(l), in terms of TT, TB, BT, and BB-states. 
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FIG. 5. Scheme to visualize the inductive proof for the strong collective decoherence case. TB- and BT-states of DFS„(J) 
are indicated. su(n J _ 1 / 2 ) ac ts on DFS„(J — 1) and DFS„(J) via DFS n -i(J — 1/2). See text for details. 
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